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Abstract 


The  first  phase  of  the  study  of  the  performance  of  a  wormgear 

transmission  is  reported.  In  this  phase  the  work  included  the  selection  of  a 

double -enveloping  wormgear  type,  and  its  dimensions,  suitable  for  use  in 
helicopter  transmissions;  the  3-D  graphics  representation  of  the  selected 
wormgear  using  the  I-DEAS  software;  the  analysis  of  the  kinematics  of  meshing; 
the  analysis  of  load  sharing  among  the  meshing  teeth;  and  the  implementation 
of  the  analyses  in  a  computer  program.  The  report  describes  the  analyses, 
their  results,  and  the  use  of  the  computer  programs. 
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1  BackgroTind  and  Objectives  of  the  Study 


The  wormgear  is  attractive  for  use  in  helicopter  transmissions  because 
it  can  accomplish  two  functions  in  one  stage,  viz.  transmitting  torque  from 
the  horizontal  engine  shaft  to  the  vertical  rotor  shaft,  and  simultaneously 
realizing  the  required  speed  ratio.  Thus,  the  use  of  the  wormgear  can  result 
in  compact  transmission  design  with  significant  weight  saving.  However,  the 
wormgear  has  a  serious  drawback.  Because  worm  and  gear  teeth  are  in  pure 
sliding  contact,  the  sliding  friction  causes  excessive  power  loss  and  renders 
the  efficiency  of  power  transmission  low.  Recently,  Chaiko  (1990)  advanced 
the  concept  of  using  hydrostatic  lubrication  to  remove  the  sliding  friction. 

He  considered  a  wormgear  transmission  of  realistic  dimensions  and  under 
realistic  operating  condition,  and  found  that  the  efficiency  of  the 
transmission  could  reach  99.14%,  and  its  weight  could  be  only  41%  that  of  the 
baseline  transmission.  These  figures  convincingly  demonstrate  that  the 
hydrostatically  lubricated  wormgear  can  be  a  candidate  of  the  advanced 
helicopter  transmissions.  To  further  investigate  the  performance  of  this  type 
of  transmission,  an  analytical  tool  capable  of  generating  quantitative  data 
for  a  wide  range  of  operating  conditions  is  needed.  Hence  a  research  proposal 
was  submitted,  and  a  grant  was  subsequently  awarded  (NAG3-1316),  to  develop 
such  a  tool. 

The  project  consists  of  three  years  of  study  with  the  following 
obj  ectives : 

1.  Year  1  is  to  be  devoted  to  the  selection  of  a  double -enveloping  wormgear 
type  suitable  for  use  in  helicopter  transmissions,  the  analysis  of  the 
kinematics  of  meshing,  and  the  analysis  of  load  sharing  among  the 
meshing  teeth. 

2.  Year  2  is  to  be  devoted  to  the  selection  of  a  hydrostatic  lubrication 
configuration  on  the  gear  surface,  and  a  comprehensive  analysis  of  the 
hybrid  (hydrostatic  and  hydrodynamic)  lubrication  at  the  contact  based 
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on  the  surface  geometry  and  running  conditions  obtained  in  Year  1 . 

3.  Year  3  is  to  be  devoted  to  the  analyses  of  accessory  bearings  in  the 

wormgear  transmission,  and  the  assembly  of  the  obtained  information  in  a 
computer  program  package  for  the  design  of  this  type  of  transmission. 

The  present  report  describes  the  work  done  in  Year  1 . 


2  Wormgear  Geometry 

2.1  Selection  of  Wormgear  Type 

Wormgears  are  used  to  transmit  power  between  two  non-parallel  non¬ 
intersecting  shafts.  The  angle  between  the  shafts  is  usually  90  deg.  There 
are  two  common  types  of  wormgears:  single -enveloping  and  double -enveloping. 
Fig.  2.1  shows  a  single-enveloping  wormgear  where  the  gear  envelops  the  worm. 
The  double -enveloping  wormgear  is  shown  in  Fig.  2.2,  which  is  like  the  single- 
enveloping  wormgear  except  that  the  worm  also  envelops  the  gear. 

The  double -enveloping  wormgear  has  two  contact  lines  on  the  mating 
surface  during  most  of  the  meshing  period.  In  comparison,  the  single- 
enveloping  wormgear  has  only  one  contact  line.  Hence,  the  double -enveloping 
wormgear  is  expected  to  have  superior  contact  properties .  The  usual  way  to 
manufacture  a  wormgear  is  first  milling  a  worm,  then  cutting  the  gear  with  a 
gear  hob  which  has  the  same  shape  as  the  worm.  The  worm  and  gear  hob  are 
usually  ground  for  high  precision  applications.  The  grinding  wheel  can  be  a 
disc  (Shen,  et  al.,  1983),  Fig.  2.3,  or  a  cone  (Simon,  1989),  Fig.  2.4.  The 
worm  grinding  process  is  called  the  first  enveloping  process;  and  the  gear 
cutting  process  is  called  the  second  enveloping  process. 

In  the  first  enveloping  process  the  grinding  wheel  actually  rotates 
around  a  base  circle,  with  angular  velocity  ^2>  wtiile  the  worm  rotates  about 
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its  own  axis  with  angular  velocity  The  distance  between  the  worm  axis  and 

»  * 

the  base  circle  axis  is  denoted  as  Sq,  and  the  ratio  ^>1  to  V2  is  denoted  as  i(j. 
In  the  case  of  the  disc  grinding  wheel,  the  grinding  plane  may  be  parallel  to 
the  base  circle  axis  or  inclined  to  it  at  an  angle  jS.  The  center  distance  a 
and  speed  ratio  i  of  the  wormgear  are  actually  the  respective  quantities  used 
in  the  second  enveloping  process.  If  a=aQ  and  i=io,  the  wormgear  is  said  to 
be  of  the  conventional  form;  otherwise  it  is  said  to  be  of  the  modified  form. 

In  helicopter  applications  the  torque  load  and  the  speed  ratio  are  high. 
To  accommodate  such  stringent  operating  conditions  the  double-enveloping 
wormgear  was  selected.  This  is  not  only  because  the  double-enveloping 
wormgear  has  a  larger  contact  region  but  also  because  its  contact  pattern  is 
favorable  to  the  establishment  of  fluid  film  lubrication.  Further,  for  ease 
of  manufacturing  the  disc  grinding  method  was  chosen.  The  analysis  presented 
below  is  based  on  the  conventional  form,  but  it  can  be  easily  extended  to  the 
modified  form. 

2.2  Layout  of  Wormgear  Dimensions 

Wormgear  dimensions  are  usually  determined  by  following  the  gear 
standards  (e.g.  AGMA,  1965),  which  were  developed  based  on  many  years  of 
experience.  The  procedure  typically  starts  with  a  known  speed  ratio  and  a 
known  center  distance.  For  helicopter  applications,  the  suggested 
transmission  power  output  is  4000  kW  ,  the  speed  ratio  is  110,  and  the  center 
distance  is  about  1500  mm.(Chaiko,  1990)  These  parameters  are  beyond  the 
range  covered  by  the  AGMA  standards.  Hence,  while  the  standards  were  followed 
as  closely  as  possible,  some  modifications  had  to  be  made.  Table  1  lists  the 
given  and  the  calculated  dimensions  and  parameters,  along  with  the  formulas 
used  for  the  calculation. 

Common  terminologies  used  to  describe  the  double-enveloping  wormgear  are 


3 


shown  in  Fig.  2.5.  Another  term  that  should  be  introduced  is  the  mid-plane, 
which  refers  to  the  plane  passing  through  the  worm  axis  and  perpendicular  to 
the  gear  axis.  The  pitch  circle  of  the  gear  is  measured  in  the  mid-plane. 

The  pitch  circle  of  the  worm  is  measured  at  the  center  of  the  worm.  In  the 
double -enveloping  wormgear  the  worm  circular  pitch  and  the  gear  circular  pitch 
are  equal.  The  base  circle  is  used  to  define  the  worm  profile  in  the  mid¬ 
plane.  The  normal  pressure  angle  is  always  chosen  to  be  20  deg. 

Several  of  the  selected  quantities  in  Table  1  should  be  explained.  A 
right  hand  worm  was  selected.  The  pitch  diameter  of  the  worm  dj^  was  chosen  to 
be  300  mm  because  of  strength  and  rigidity  considerations.  The  number  of 
meshing  teeth  Np  was  chosen  to  be  10  instead  of  the  usual  4  or  5  in  order  to 
reduce  the  load  level  on  the  meshing  teeth.  The  choice  of  ^  =  5  deg  was  a 
compromise  between  acquiring  a  better  contact  condition  and  maintaining  an 
acceptable  tooth  thickness  at  the  top  of  the  worm  thread. 

It  should  be  noted  that  for  the  double -enveloping  wormgear  a  different 
choice  of  dimensions  can  result  in  very  different  surface  shape  and  contact 
pattern.  The  results  shown  below  are  based  on  the  data  of  Table  1. 

To  better  visualize  the  selected  wormgear,  a  software  named  I-DEIAS  was 
used  to  generate  its  3-D  graphics.  Figs.  2.6  to  2.9.  Fig.  2.6  shows  the 
engaged  wormgear.  The  worm  and  gear  are  shown  separately  in  Figs.  2.7  and 
2.8.  Fig.  2.9  depicts  a  close-up  view  of  the  gear  teeth. 

Details  of  the  I-DEAS  program  that  generated  these  graphics,  and  the 
procedure  to  run  the  program,  will  be  described  in  Section  5.1. 
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3  Contact  analysis 


3.1  Method  of  Enveloping 

The  geometric  and  kinematic  properties  of  the  wormgear  were  analyzed 
with  the  method  of  enveloping.  To  describe  the  method,  let  us  consider  two 
rotating  machine  components,  m  and  k,  which  are  in  meshing.  Let  us  attach  a 
coordinate  system  S„(x„,y„,z„)  to  the  machine  component  m;  and  likewise  attach 
a  coordinate  system  Sk(Xk,yk,Zi.)  to  the  machine  component  k.  The  two  machine 
components  are  rotating  with  the  angular  coordinates  and  Vk  with  respect  to 
some  stationary  reference  frame.  Since  the  two  components  are  in  meshing, 
there  must  be  some  definite  relationship  between  V’m  Hence  only  one  of 

them  is  independent,  and  either  can  serve  as  the  running  parameter.  Let  V 
denote  this  running  parameter.  The  surface  of  machine  component  m  can  be 
expressed  as  f (x„,ym,z^)=0  in  S„.  By  coordinate  transformation,  f=0  becomes  a 
family  of  surfaces,  F(Xk,y]j,Zi,,V')=0  in  Sj,.  Then  the  surface  of  the  machine 
component  k  is  the  envelope  of  F^O,  viz. 

(1) 

3F/3i|r=0 

If  =  constant,  the  surfaces  of  the  two  meshing  components  are  said  to  be 

conjugate  surfaces  (Litvin,  1989) . 

While  the  idea  was  simple,  the  algebraic  manipulations  needed  to  carry 
it  out  were  quite  clumsy.  Recently,  Shen,  et  al.  (1983)  reformulated  the  idea 
using  a  set  of  operators.  Because  many  of  the  operators  involving  the 
transformation  of  coordinates  can  be  done  once  for  all,  the  method  has  become 
simple  as  well  as  elegant.  In  the  special  case  of  conjugate  surfaces,  the 
method  is  even  more  convenient  to  use.  This  will  be  demonstrated  below.  The 
reader  is  referred  to  Shen,  et  al.  (1983)  for  the  details  of  the  method. 
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3.2  Coordinate  Systems  and  Operators 


To  carry  out  the  analysis,  four  right  hand  coordinate  systems,  Sj^,  S2, 

Sj  and  Sjj,  are  set  up.  See  Fig.  3.1. 

Sj  and  Sj,  are  stationary  systems.  Sj(Xj,yj,Zj)  is  located  at  the  worm 
position.  Xj  is  along  the  center  line  of  the  wormgear  and  points  to  the  gear, 
and  Zj  axis  is  along  the  worm  axis.  S„(x^,yjj, z„)  is  located  at  the  gear 
position.  Xjj  is  along  the  center  line  but  opposite  to  Xj ,  and  z^  is  along  the 
gear  axis . 

S],  and  S2  are  moving  systems.  Si(Xi,yj^,Zi)  rotates  with  the  worm  and  its 
z^  axis  always  coincides  with  the  Zj  axis.  The  worm  rotation  angle  is  the 
angle  between  Xj^  and  Xj .  S2(x2,y2,Z2)  rotates  with  the  gear  and  its  Z2  axis 

always  coincides  with  the  axis.  (In  the  first  enveloping  process  S2  turns 
around  the  base  circle  with  the  grinding  wheel.)  The  rotation  angle  ^^2 
angle  between  X2  and  x^.  Thus,  is  fixed  to  the  worm  and  S2  is  fixed  to  the 
gear;  and  the  worm  and  gear  rotate  about  Sj  and  Sjj,  respectively.  Initially 
Si  coincides  with  Sj ,  and  S2  coincides  with  Sj,. 

Let  be  the  coordinate  transformation  matrix  from  Sj„  to  Sj.,  where 

the  subscript  (  indicates  that  V  is  the  running  parameter  of  the  coordinate 
transformation.  If  the  running  parameter  is  understood,  the  subscript  is 
dropped  to  simplify  the  notation.  In  addition,  if  the  transformation  is 
between  two  stationary  systems,  then  no  running  parameter  is  involved  and 
there  will  be  no  subscript.  The  coordinates  in  can  be  transformed  into  Sj. 
as : 


Yk 

=  a;*”’ 

Ym 

Zk 

1 

1. 

(2) 
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and  a  function  in  S„  can  be  transformed  into  Sj.  as: 

(x„.y„,z„)  =f  (3) 

The  coordinate  transformation  matrices  from  Sj  to  Sj,  from  S2  to  S„,  and 
from  Sj,  to  Sj  are  given  as: 

cosilfi  simjri  0  O' 

Aiij)  =  -sin’I'i  cosTir,  0  0  (4) 

0  0  10 
0  0  0  1 


cosilTj  -sim(r2  0  O' 

sintjra  costlr^  0  0  (5^ 

0  0  10 
0  0  0  1. 


-10  0  do 

0  0-10  (6) 
0-100 
0  0  0  1, 

The  transformation  matrix  from  S2  to  may  be  written  as 
or: 

-cosi|rj.cosi|r2  cosi|;iSini|r2  -sini|ri  ao*cosi|ri 
sinijriCosilrj  -sinilJiSini|»2  -cosi|ri  -ao^simlTj 
-sinilTj  -cos^lfj  0  0 

0  0  0  1 

Note  that  the  relation  a'‘*’=(A^*°“^)'^  holds. 


In  addition  to  the  coordinate  transformation  matrix,  two  other  operators 
are  defined.  One  is  the  differentiation  operator: 
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The  other  is  the  similar  differentiation  operator: 


In  the  above  formulas,  the  subscript  (  may  be  dropped  if  no  confusion  can 
arise . 


Several  needed  properties  of  these  operators  are  given  below; 

=A<*">Df+  PA<*">y„+ 


i2AW)f=;!,W>j9Wf 


(11) 


(12) 


where 


(13) 


In  the  above  formulas,  k  and  m  can  be  1,  2,  j  and  n;  but  i  and  j  can  only  be  1 
and  2 . 


To  distinguish  between  the  first  and  the  second  enveloping  processes,  ^ 
is  used  as  the  running  parameter  for  the  first  enveloping  process  and  S  is 
used  for  the  second  enveloping  process.  The  order  of  operations  with  respect 
to  different  running  parameters  may  be  interchanged,  e.g. 


(14) 


The  derivation  of  the  above  and  other  properties  of  the  operators  can  be  found 
in  Shen,  et  al.  (1983). 

It  should  be  noted  that  different  symbols,  Iq  and  i,  are  used  to  denote 
the  respective  speed  ratios  in  the  first  and  the  second  enveloping  processes, 
i .  e . 


h 

^2 

h 

02 


(15) 


As  said  before,  if  ig  =  i  and  ao=a  the  wormgear  is  of  the  conventional  form; 
otherwise  it  is  of  the  modified  form. 


3 . 3  Worm  Surface 

The  worm  surface  is  formed  by  a  grinding  wheel  in  the  first  enveloping 
process.  The  surface  of  the  grinding  wheel  is  called  the  original  generating 
surface.  For  a  disc  grinding  wheel  this  original  generating  surface  is  a 
plane.  For  a  cone  grinding  wheel  the  original  generating  surface  is  a  line  on 
the  cone. 


The  original  generating  surface  can  be  expressed  in  the  S2  system  as 
f (x2,y2,Z2)=0.  In  the  first  enveloping  process  S2  rotates  about  Z2  with  the 
rotation  angle  ^2.  rotates  about  Zj  with  the  rotation  angle  Vi-  As  a 

result  of  rotation,  the  original  generating  surface  forms  a  one -parameter 
family  in  S^,  F(xi,yi,Zi,^)=A^‘^2’f (x2,y2.Z2)=0.  Because  ^1/^2  =  constant,  the 
running  parameter  ip  can  be  either  V*!  ^2'  Then,  the  envelope  of  the  family 


9 


of  surfaces  is  the  worm  surface,  which  is  given  as  F=0  and  D^F=0. 

Considering  only  the  case  of  the  selected  wormgear,  the  original 

generating  surface  is  a  plane  in  S2,  which  is: 

f  {x2,y2,Z2)=y2-Z2tanP+ri,=  0  (16) 

After  coordinate  transformation  into  Sj^,  the  one-parameter  family  of  planes 
becomes : 

f=x.^  (cosilfj^sinijrj+tanPsimlfi) 

+yi  (-sinf  j^sin\|;2+tanPcosi|ri) 

-Z2cosij;2-(aosini|r2-rj,)  =0 


Then,  the  worm  surface  may  be  represented  by  the  following  two  equations: 

F  =  Xj  (cositfisinijfj+tanpsinilfi)  +  yj  (-sini|riSin»|r2+tanPcosT|r2) 
-2iCosilr2-(aoSini(r2--ti)  =  0 
D^F  =  jf^CcosilTiCOSTlTj+iotanpcosilfi-ioSinifiSinTlrj) 

+  y^  (-simjr^cosilrj-iotanpsinti-iocosifisinilrj)  +Z2sini|r2-a(,cost|r2  =  0 


Using  the  relation  eliminating  Vi  and  V’2  from  Eqs.  18,  an  explicit 

equation,  which  is  nonlinear,  of  the  worm  surface  can  be  obtained. 

For  a  given  value  of  Eqs.  18  represent  the  equation  of  a  straight 
line.  Recall  that  the  first  equation  is  the  equation  of  the  original 
generating  plane  expressed  at  the  position  Hence  the  straight  line  is  on 

both  the  generating  plane  and  the  enveloping  surface.  This  straight  line  is 
the  contact  line  between  the  worm  and  the  disc  grinding  wheel  at  the  position 
^>1.  It  can  also  be  seen  that  the  worm  surface  is  the  set  of  all  the  straight 
contact  lines.  Fig.  3.2  depicts  this  situation. 
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According  to  differential  geometry,  if  the  planes  in  a  one -parameter 
family  are  not  parallel  and  there  is  not  a  point  common  to  all  the  planes, 
then  the  envelope  of  this  family  is  the  tangent  surface  of  its  edge  curve 
(Kreyszig,  1959).  In  Fig.  3.2  the  curve  labelled  T  is  the  edge  curve.  The 
edge  curve  is  given  by  these  equations: 


F  =  Xi(cosi|riSinilr2+tanPsini|ri+yi(-simIriSini|r2+tanpcosi|ii-ZiCOST|r2-(ao*sini|r2-r^)  =  0 
D,F=Xi(cosi|fiCOS»|r2+iotanPcosilri-ioSini|f2sin»iJ2) 

+  yi(-sin^ricosilr2-iotanpsin\|ri-i(,cosi|risini|r2)  +ZiSini|r2-aoCosilr2=0 
D^D^F=  (iocos»lrisini|i2+iotanpsini|ji+2*iosini|ricosi|f2+cosijrisim|r2)  ^3^9) 

+  yi(-ioSini|riSinilr2+iotanPcosT|ri+2*ioCosilriCosi|>2-sini|riSini|r2)  -ZiCOSTlr2-aoSini|r2=  0 
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The  edge  curve  divides  the  envelope  surface  into  two  branches,  of  which 
one  is  the  worm  surface.  If  the  worm  thread  domain  determined  by  the  selected 
dimensions  intersects  the  edge  curve,  the  surface  of  the  disc  grinding  wheel 
would  extend  over  the  edge  curve.  The  crossed- over  portion  of  the  grinding 
surface  would  undercut  the  worm  thread.  Undercutting  the  thread  would 
decrease  the  length  of  contact  lines,  decrease  the  root  thickness  of  the  worm 
thread,  and  consequently  decrease  the  load  carrying  capacity  of  the  thread. 
Hence,  wormgear  dimensions  should  be  chosen  properly  so  that  the  worm  thread 
domain  is  away  from  the  edge  curve. 

3.4  Gear  Surface 

The  same  procedure  applies  to  the  second  enveloping  process.  In  this 
process  the  gear  hob  has  the  same  shape  as  the  worm  surface.  Hence,  the 
generating  surface  is  given  by  Eqs.  18  in  S^.  As  the  result  of  rotation,  this 
generating  surface  forms  a  one -parameter  family  of  surfaces  in  S2.  The 
running  parameter  for  this  process  can  be  either  or  62,  while  becomes  a 
geometric  parameter  already  in  existence  in  the  worm  surface  equation.  The 
procedure  to  obtain  the  gear  surface  may  be  as  follows:  By  eliminating  ip  from 
Eqs.  18,  one  obtains  an  explicit  equation  of  the  worm  surface.  This  equation 
may  be  written  as ; 

F*(Xi,yi,Zj)  =0  (20) 

The  one -parameter  family  of  surfaces  in  S2  is  obtained  by  coordinate 
transformation : 

(X2,y2,zj,e)=  F*=  0 

The  envelope  of  the  family: 

f  <2)=  0 

Df  F  =  0 


(21) 


(22) 


is  then  the  gear  surface. 
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But  it  is  not  easy  to  obtain  F*=0.  Neither  is  it  necessary.  Assvime 
that  V^V’(Xi,yi,Zi)  is  obtained  by  solving  the  second  equation  of  Eqs.  18, 

D^F=0.  Substituting  it  into  the  first  equation  of  Eqs.  18,  F=0,  yields  the 
explicit  worm  surface  equation 

F*=  F  0  (23) 

Then,  the  equation  Dgf can  be  written  as: 

(24) 

But  from  Eq.  23, 

<25) 

Then,  since  D^F=0  and  Eq.  23  does  not  contain  6,  Eq.  24  reduces  to: 

which,  on  multiplification  by  Ag'^^^Ag^^^\  becomes  Ae^^^^Dg^^^F=0 .  Thus,  the  gear 


surface  is  given  by: 


=  0 

^(21)p^p.  ^  0 

j^i2DDii>F  =  0 


which  can  be  written  as: 


*11  ®12  “13  “14 


®21  ^22  ^23  ^24  Fa  _ 


where 
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(29) 


aii=-sinilr2cos02cos  (+1-61)  -tanPcosB^sin  +cosijr2sine2 

33^2“  sini|r2sin02cos  {i|ri-02)  +tanpsin02sin{T|r3-0i)  +cos\|r2cos02 
333=  sinilr2sin(ilr3-03)  -tanPcos  (ijr^-S^) 

3^4=  a  [simlTjCOS  (itri-6i)  +tanpsin(i(r3^-03)  ]  -aoSinijrj+r^ 


323=  sini|r2cos02sin  (i(f3-03)  -tanPcos02Cos  (+3-03) 

-  [sini|r2sinO2+cosi|f2Cos02Cos  (tjfj-Oi)  ]  /io 
a22=-sinijf2sin02sin  (1^3-03)  +tanPsin02COS  (i|r3-03) 

-  [sini|r2cos02-cosilr2sin02cos  (itf3-03)  ]  / io 
323=  sinijrjcos  (1^3-03)  +tanpsin(i|f3-03)  +  [cosi^jSin  (i(r3-03)  ]  / 

324=-a  [sirnjf2sin(»|r3-03)  -tanpcos  (\|r3-03)  -cosT|r2Cos  (1^3-03)  /io]  -aocostlf2/io 


333=  a32-i*  [sinilf2cos02sin (11^3-03)  -tanPcos02Cos  (1^3-03)  ] 
332=-a33+i*  [sini|r2sin02sin(ijr3-03)  -tanPsinOjCos  (!ji3-03)  ] 
a33=-i*  [sinipjcos  (»|J3-03)  +tanPsin (4r3-03)  ] 

334=  a*i*  [sinT|r2sin(i|r3-03)  -tanPcos  (+3-03)  1 


In  Eqs.  27  there  are  three  equations  and  five  unknowns.  Using  the 
relations  Vi/V’2  =  ^1/^2  explicit  form  of  the  envelope  can  be 

obtained  by  eliminating  tj)  and  0.  Eqs.  27,  along  with  the  relation  limiting 
the  gear  tooth  domain,  then  provide  the  shape  of  the  gear  tooth. 

Consider  the  case  of  the  modified  form.  For  a  given  value  of  6,  Eqs.  27 
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produce  a  curve.  The  curve  is  the  contact  line  between  the  worm  and  the  gear 
at  this  ff  position.  The  curve  goes  over  the  gear  surface  twice  as  shown  in 
Fig.  3.3.  For  fixed  values  of  ip  and  6,  Eqs.  27  represent  a  point  on  the 
contact  line. 

Consider  the  case  of  the  conventional  form.  For  a  given  6 ,  if  one 
chooses  ■>p=6 ,  the  third  and  the  second  equations  of  Eqs.  27  are  identical. 

Then  Eqs.  27  reproduce  a  straight  contact  line  of  the  first  enveloping 
process.  This  kind  of  contact  lines  is  called  the  first  enveloping  contact 
lines  (the  solid  lines  in  Fig.  3.4).  The  set  of  all  these  straight  lines 
reproduce  a  part  of  the  original  generating  plane.  For  a  given  6,  let 
Eqs.  27  give  a  curved  contact  line,  called  the  second  enveloping  contact  line 
(the  dashed  lines  in  Fig.  3.4).  The  set  of  all  the  curved  contact  lines  form 
the  curved  part  of  the  gear  surface .  Thus ,  in  the  case  of  the  conventional 
form,  the  gear  surface  consists  of  a  planar  region  and  a  curved  region.  For  a 
certain  range  of  6  values  these  two  regions  overlap  each  other.  The 
overlapped  region  is  shown  in  Fig.  3.4  where  the  dashed  lines  and  the  solid 
lines  cross.  This  situation  is  impossible  to  occur  in  reality.  From  a  view 
of  the  (72,22)  plane  it  is  clear  that  the  curved  region  is  farther  away  from 
the  solid  interior  than  the  planar  region.  Hence,  the  curved  region  is 
fictitious,  which  is  actually  bobbed  off  during  manufacturing.  Analytically 
the  principal  relative  curvature  is  used  to  sort  out  the  fictitious  part 
(Shen,  et  al,  1983).  The  relative  curvature  is  the  difference  of  curvatures 
of  the  worm  and  gear  surfaces  at  the  contact  point.  Because  the  gear  tooth 
envelops  the  worm  thread,  at  the  contact  point  the  curvature  of  the  gear  tooth 
must  not  be  greater  than  that  of  the  worm  thread  (absolute  values).  Hence,  if 
the  relative  curvature  is  negative ,  the  produced  gear  surface  is  not  real . 

The  derivation  of  the  expression  for  the  principal  relative  cuirvature  is  given 
in  Appendix  3. 

From  the  above  description  it  can  be  seen  that,  in  the  double -enveloping 
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wormgear,  the  worm  and  gear  are  in  line  contact  and  there  are  usually  two 
contact  lines  on  each  pair  of  the  meshing  teeth. 

3.5  Tooth  domains 

After  the  equations  of  the  worm  and  gear  surfaces  are  obtained,  the 
tooth  domains  need  be  specified  to  determine  what  portions  of  these  surfaces 
are  the  actual  worm  thread  and  gear  tooth.  Table  1  can  be  used  to  set  up  the 
following  geometrical  limitations  on  the  worm  and  gear  surfaces.  The  worm  and 
gear  tooth  domains  are  shown  in  Fig.  3.5.  Several  new  symbols  that  appear  in 
the  following  equations  are  also  defined  in  this  figure. 


9 

where  iCy,  is  the  radius  of  any  point  on  the  worm  surface . 


The  maximum  and  minimum  radii  of  the  gear  tooth  are; 
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R^=a-^{dj2V-(b^/2)- 


r^=  min  {(a+Aa)  -  )fSi^-Tz^/2)^ ,  } 

rg^=(a+Aa)  - 

r,=  (33) 

^^5-x 


where  rg  is  the  radius  of  any  point  on  the  gear  surface,  and 


a  -  ^{d,j2V-^bJ2V 

4,i„=  a  -  ^<,dj2V-{bJ2)-  (34) 

Aa  =  {h^*dfj2)  -  a 


3.6  Properties  Needed  for  Lubrication  Studies 

To  carry  out  the  lubrication  analysis  several  geometric  and  kinematic 
pj;op©i’ties  are  needed.  These  include  the  shape  of  the  two  contacting 
surfaces,  the  relative  velocity  between  the  two  surfaces  at  a  contact  point, 
and  the  normal  forces  acting  on  the  meshing  teeth,  etc. 

3.6.1  Surfaces  in  Stationary  Systems 

Eqs.  18  and  Eqs.  27  are  obtained  in  the  moving  coordinate  systems,  where 
the  worm  and  the  gear  appear  stationary.  To  see  the  moving  worm  and  gear 
surfaces ,  these  equations  must  be  expressed  in  the  stationary  systems .  At  a 
given  instance,  or  for  a  given  gear  position  9,  there  are  two  ways  to  obtain 
the  worm  and  gear  surfaces  in  a  stationary  coordinate  system.  One  way  is 
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calculating  the  coordinates  in  the  moving  system  first ,  then  transforming 
these  points  into  the  stationary  system.  The  other  way  is  transforming  the 
equations  into  the  stationary  system  first.  In  either  way  the  coordinate 
transformation  matrices  are  specialized  to  the  given  position.  To  distinguish 
this  position  from  the  running  parameter  6,  it  is  denoted  as  9*.  Thus  the 
worm  surface  in  Sj  is  given  as : 


=  0 
=  0 


(35) 


the  gear  surface  in  Sj^  is  given  as: 


=  0 

A^^^'^A^^^'^D^F  =  0 

^<n2)^(2i)^(i)^  =  0 


and  the  gear  surface  in  Sj  is  given  as: 


(37) 


3.6.2  The  Relative  Velocity  at  a  Contact  Point 

In  Eqs.  35  to  37  if  one  sets  9=9*,  the  contact  lines  expressed  in  the 
stationary  systems  are  obtained  for  the  given  9*  position.  Then  the  relative 
velocity  at  a  contact  point  can  be  calculated,  which  is  the  difference  of  the 
worm  velocity  and  the  gear  velocity  at  this  contact  point.  In  Sj  the  relative 
velocity  is  given  as: 
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(38) 


Vl=-yj  T  +  Xj  j 

^=-yr“2<-i  )  +  ) 

where  Xj,  y^,  Zj  are  the  coordinates  of  the  contact  point  in  Sjj  Xj,.  yn> 
the  corresponding  coordinates  in  S^;  i,  j  and  k  are  the  unit  vectors  in  Sj . 

To  facilitate  further  development,  let  us  introduce  the  term  "center 
point" .  A  center  point  is  the  point  of  intersection  between  a  tooth  surface 
and  the  gear  pitch  circle.  In  the  meshing  zone  there  is  a  gear  center  point 
on  every  gear  tooth.  Correspondingly  there  is  a  worm  center  point  on  the  worm 
coil  in  meshing  with  this  gear  tooth.  The  worm  center  point  is  in  general 
slightly  different  than  the  gear  center  point,  unless  the  center  point  is  a 
contact  point.  In  the  analysis  below,  the  normal  and  friction  forces  between 
a  pair  of  meshing  teeth  are  taken  as  acting  at  the  worm  center  point,  which 
for  brevity  will  be  referred  to  simply  as  the  center  point. 

In  case  that  the  center  point  is  a  contact  point,  the  relative  velocity 
is : 


Vr  =  =  y^Uj.Ci>i)2+(c4«2/2)* 

yj  =  0,  +  yn  =  <  )* 


3.6.3  Normal  Direction  of  the  Worm  Surface 

Because  the  worm  surface  is  formed  by  the  original  generating  plane,  the 
normal  direction  at  a  point  of  the  worm  surface  is  the  same  as  the  normal 
direction  of  the  original  generating  plane.  Further,  since  the  worm  surface 
and  the  original  generating  plane  contact  in  a  straight  line,  all  the  points 
on  the  contact  line  have  the  same  normal  direction.  In  this  normal 
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direction  is: 


njf=  cosi|rjSinT|r2+tanPsini|ri 
xiy=-siniJriSini|r2+tanPcosT|rj^ 

JI^=-C0St2 


(40) 


Transforming  Eqs .  40  into  $2,  the  normal  direction  of  the  gear  surface 
at  the  contact  point  can  also  be  obtained.  Note  that  since  the  normal 
direction  of  the  worm  is  outward  from  the  worm  body,  the  normal  direction  of 
the  gear  is  then  inward  to  the  gear  body.  The  normal  direction  can  be  further 
transformed  into  a  stationary  coordinate  system.  Thus,  for  the  given  position 
6* ,  the  normal  direction  in  Sj  is: 

cos  (ijrj-6i)sini|r2+  tanPsin 

njy=-sin  (il»2-0i)  sin»|J2+  tanPcos  (41) 

iijz=-cosi|r2 


3.6.4  Clearance  Function 

Once  the  worm  and  gear  surfaces  are  found,  the  clearance  between  the  two 
mating  surfaces  can  be  determined.  This  information  is  of  vital  importance  in 
the  analysis  of  fluid  film  lubrication.  This  clearance  is  defined  to  be  the 
gap  between  the  two  surfaces  measured  in  the  direction  of  the  common  normal  on 
the  contact  line.  Clearly,  as  the  distance  from  the  contact  line  increases, 
the  gap  so  measured  becomes  less  representative  of  the  true  distance  between 
the  two  surfaces.  But  this  is  immaterial,  because  as  the  gap  becomes  large, 
the  effects  of  fluid  film  lubrication  also  become  negligible.  The  clearance 
function  can  be  obtained  numerically  as  follows:  (1)  Choose  the  grid  system 
on  the  gear  tooth  and  obtain  the  coordinates  points. 

(2)  Calculate  the  normal  direction  at  a  contact  point  and  obtain  the  equation 
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of  the  straight  line  passing  through  a  grid  point  and  in  this  normal 
direction.  (3)  Find  the  point  of  intersection  between  this  straight  line  and 
the  worm  surface .  The  gap  at  the  grid  point  is  then  the  distance  between  the 
grid  point  and  the  point  of  intersection.  It  is  convenient  to  carry  out  this 
calculation  in  the  lubrication  analysis,  where  the  grid  system  has  to  be 
selected  first. 

3.7  Results  and  Discussion 

Sample  results  were  computed  for  the  selected  wormgear.  Fig.  3.6  shows 
the  wormgear  position  at  6  =0,  where  ten  pairs  of  teeth  are  in  meshing, 
beginning  from  gear  tooth  No. 2  and  ending  with  gear  tooth  No. 11.  The  contact 
lines  on  several  of  these  teeth  are  shown  in  Figs.  3.7  as  viewed  in  the  (Xj,^ 
Zjj)  plane.  In  each  of  these  figures  the  circles  are  on  the  gear  tooth,  whose 
domain  is  shown  by  the  solid  lines ;  the  stars  are  on  that  portion  of  the  worm 
thread  that  overlapps  the  gear  tooth.  Note  also  that  the  unit  of  the 
coordinates  in  these  figures  is  millimeter,  and  the  teeth  are  shown  in  their 
true  proportion.  On  gear  tooth  No. 2,  Fig.  3.7(a),  the  two  contact  lines  are 
far  from  the  mid-plane  (Zn=0) .  These  contact  lines  move  toward  the  mid-plane 
as  the  tooth  ntunber  increases.  See  Figs.  3.7(a)  to  3.7(d).  On  gear  tooth 
No. 10,  Fig.  3.7(e),  the  second  enveloping  contact  line  crosses  over  the  first 
enveloping  contact  line.  The  crossed-over  part  is  nonexistent,  as  is 
explained  previously.  On  gear  tooth  No. 11,  Fig.  3.7(f),  the  entire  second 
enveloping  contact  line  falls  below  the  first  enveloping  contact  line,  and 
becomes  non-existent.  Thus,  toward  the  end  of  the  meshing  zone,  the  total 
length  of  the  contact  lines  diminishes.  From  these  figures  it  can  also  be 
seen  that  there  are  two  contact  lines  on  the  majority  of  the  teeth  in  the 
meshing  zone;  these  contact  lines  are  located  rather  symmetrically  on  the  two 
sides  of  the  mid-plane  and  they  run  across  the  tooth  face  more  or  less  in  the 
direction  of  the  X2  axis . 
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The  total  length  of  the  contact  lines  on  each  meshing  tooth  is  plotted 
in  Fig.  3.8.  The  situation  corresponding  to  the  position  6*’=0  is  shown  in 
Fig.  3.8(a).  When  the  gear  rotates  one  degree,  i.e.  5*=1  deg,  tooth  No. 11 
moves  out  of,  while  tooth  No.l  has  not  moved  into,  the  meshing  zone.  Hence 
there  are  only  nine  pairs  of  teeth  in  meshing.  The  distribution  of  the 
contact  line  lengths  in  this  situation  is  shown  in  Fig.  3.8(b).  When  the  gear 
rotates  one  more  degree,  i.e.  d*=2  deg,  tooth  No.l  has  moved  into,  while  tooth 
No. 10  still  remains  in,  the  meshing  zone.  There  are  again  ten  pairs  of  teeth 
in  meshing.  The  corresponding  distribution  of  the  contact  line  lengths  is 
shown  in  Fig.  3.8(c).  Since  the  angle  between  two  neighboring  teeth  is  36O/N2 
=3.273  deg,  if  the  gear  further  rotates  1.27  deg,  the  situation  of  Fig.  3.8(a) 
re -appears . 

The  variation  of  the  contact  line  length  on  one  gear  tooth  as  it  moves 
through  the  meshing  zone  is  shown  in  Fig.  3.9,  which  is  essentially  the  three 
graphs  of  Figs.  3.8  combined. 

The  radial  position  of  the  edge  curve  on  the  worm  surface  is  shown  in 
Fig.  3.10.  For  comparison,  the  radial  position  of  the  hour- glass -shaped  worm 
pitch  cylinder  is  plotted  as  the  dashed  curve.  Since  the  tooth  depth  is  about 
20  mm,  it  can  be  seen  that  the  edge  curve  is  located  far  below  the  worm 
thread.  Hence  there  is  no  problem  of  undercutting  the  worm  thread. 

From  the  above  results  one  can  identify  the  following  merits  of  the 
selected  wormgear: 

1.  On  the  majority  of  the  meshing  teeth  there  are  two  contact  lines.  This 
means  that  the  contact  region  on  a  tooth  of  the  double-enveloping 
wormgear  is  larger  than  that  of  the  single-enveloping  wormgear,  which 
only  possesses  one  contact  line  on  a  meshing  tooth. 

2.  The  contact  lines  run  across  the  tooth  face  more  or  less  in  the 
direction  of  the  X2  axis.  Since  the  relative  velocity  of  the  meshing 
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surfaces  is  nearly  perpendicular  to  the  X2  axis ,  the  two  directions  meet 
at  a  large  angle.  This  feature  is  favorable  to  the  establishment  of 
fluid  film  lubrication  at  the  contact. 

3.  On  the  majority  of  the  meshing  teeth,  the  pattern  of  the  contact  lines 
is  nearly  S3ramietrical  with  respect  to  the  center  point.  This  feature 
allows  a  major  simplification  to  be  made  in  the  load  sharing  analysis. 

4.  With  the  chosen  gear  dimensions  there  is  no  problem  of  undercutting  the 
worm  thread. 

4  Load  Sharing  Analysis 

The  methodology  of  determining  how  the  load  is  distributed  among  meshing 
teeth  is  not  well  developed.  Rational  analyses  are  rare  and  complicated.  For 
spur  and  helical  gears,  the  load  sharing  analyses  are  usually  based  on  the 
variation  of  tooth  stiffness  along  contact  lines  (Elkholy,  1989).  In  the 
practical  design  of  these  gears,  the  approach  is  to  analyze  an  idealized 
meshing  model,  then  apply  various  factors  to  modify  the  results.  To  estimate 
the  magnitudes  of  these  load  distribution  factors,  the  common  approach  is  to 
assume  that  gear  teeth  behave  like  a  series  of  thin  slices  that  deflect 
independently  (Steward,  1989).  It  is  doubtful  that  these  methods  are  suitable 
to  the  double-enveloping  wormgear,  because  its  contact  pattern  is  vastly 
different  than  those  of  the  spur  and  helical  gears.  Hence,  for  the  double - 
enveloping  wormgear,  there  does  not  exist  a  theory  of  load  sharing.  Neither 
are  there  adequate  experimental  data  to  guide  the  analysis  and  design.  In  the 
following,  we  shall  base  our  analysis  of  load  sharing  on  a  hypothesis. 

4.1  Hypothesis  on  Load  Sharing 

At  a  contact  point  the  forces  between  the  meshing  surfaces  consist  of  a 
normal  component,  that  acts  in  the  direction  of  the  common  normal,  and  a 
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friction  component.  The  latter  acts  in  the  direction  of  the  relative  velocity 
and  has  a  magnitude  equal  to  the  normal  component  multiplied  by  the  Coulomb 
friction  coefficient.  Since  the  direction  of  the  normal  and  that  of  the 
relative  velocity  vary  with  position,  these  forces  vary  from  contact  point  to 
contact  point  not  only  in  magnitude  but  also  in  direction.  However,  the 
contact  analysis  in  Section  3  reveals  that  the  contact  line  pattern  on  each 
tooth  is  nearly  symmetrical  with  respect  to  the  center  point.  In  addition, 
because  of  the  large  speed  ratio  (and  consequently  the  small  lead  angle)  of 
the  selected  wormgear,  the  variation  of  the  normal  direction  on  a  tooth 
surface  is  small.  These  two  features  allow  the  following  simplifications  to 
be  made : 

1 .  The  resultant  normal  force  on  a  tooth  is  proportional  to  the  contact 
line  length  on  the  tooth.  It  acts  at  the  center  point  and  is  in  the 
normal  direction  of  the  center  point. 

2 .  The  resultant  friction  force  on  a  tooth  acts  at  the  center  point  and  is 
in  the  direction  of  the  relative  velocity  at  the  center  point. 

The  above  simplifications  may  be  considered  as  a  hypothesis  of  load 
sharing,  whose  validity  should  be  verified  later  when  experimental  data  or 
results  of  a  deformation  analysis  become  available. 

4.2  Normal  Forces  on  Gear  Teeth 


The  coordinates  (Xn,p,z„^)  in  and  the  parameter  of  the  worm  center 
point  for  the  mth  pair  of  meshing  teeth  can  be  found  by  solving  Eqs .  35  and 
the  following  set  of  equations  describing  the  gear  pitch  circle: 


(a-Xj)^  +  zj  =  (  dj2 


(42) 


The  normal  direction  at  the  center  point  can  be  calculated  from  Eqs.  41,  where 
is  to  be  used  for  ip,  and  is  the  given  position. 
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From  Fig.  4.1,  where  friction  is  not  shown,  the  torque  contributed  by 
the  normal  force  on  the  gear  tooth  of  the  mth  meshing  pair  is: 

To- 

a-.,, 

Ilff= 

where  r„„  is  an  effective  radius;  n^^^,  n^y,  n^^  are  as  given  in  Eqs.  41;  [  ]„! 


denotes  the  mth  meshing  pair. 


By  hypothesis,  Fq^”  is  proportional  to  the  length  of  contact  lines. 
Besides,  TG=STgfc.  Hence  the  formula: 


C*"  — 


where 


F  "  is  the  normal  force  acting  on  the  gear  tooth  of  the  mth  meshing 


pair. 


Tg  is  the  required  output  torque  from  the  gear  shaft, 

is  the  length  of  contact  lines  on  the  mth  pair  of  the  meshing 
teeth,  and 

Np'  is  the  actual  niimber  of  meshing  pairs. 


The  relative  motion  between  the  worm  and  gear  is  pure  sliding.  Hence, 
friction  plays  an  important  role  in  the  performance  of  the  wormgear  (Shigley 
and  Mitchell,  1983).  Let  n  be  the  Coulomb  friction  coefficient,  then  the 
magnitude  of  friction  is  and  its  direction  on  the  gear  is  the  same  as 

the  direction  of  the  relative  velocity  V1-V2. 


At  the  center  point  of  the  mth  pair  the  worm  velocity  is  parallel  to  the 
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(Xj.Yj)  plane  and  equal  to  Xn^Wi,  the  gear  velocity  is  in  the  (xj.Zj)  plane  and 
equal  to  Sj  the  direction  of  the  relative  velocity  is  given  by 

Eqs.  39  and  is  shown  in  Fig.  4.2. 


4.3  Tangential  Forces  on  the  Worm 

The  force  acting  on  the  worm  can  be  decomposed  into  the  radial, 
tangential,  and  axial  components,  viz.  i  and  (Shigley  and 

Mitchell,  1983).  See  Fig.  4.1.  The  tangential  component  balances  the  input 
worm  torque,  which  in  turn  determines  the  input  power.  The  ratio  of  the  input 
power  without  friction  to  that  with  friction,  averaged  over  an  interval  where 
the  gear  rotates  3.273  deg,  determines  the  efficiency  of  the  wormgear. 


The  tangential  force  on  the  worm  for  the  mth  pair  of  the  meshing  teeth 
is : 

without  friction 


26 


with  friction 


F^=  -Ffoi,  [ 


l%y| 


J  (  +  (  d,©,  /2)^ 


n,. 


(48) 


Since  the  worm  torque  equals  the  sum  of  the  torques  contributed  from  the 
meshing  teeth,  it  is: 
without  friction 


ni, 


Tff  ~ 


ffl-i 


(49) 


with  friction 


Tfff  =  ^  Ff4t 
m-1 


(50) 


4.4  Results  and  Discussion 

Load  sharing  was  computed  for  a  power  output  of  4000  kW  and  a  rotation 
speed  of  130  r/min  of  the  gear  shaft.  To  account  for  the  effect  of  friction 
at  the  worm  and  gear  contact,  a  Coulomb  friction  coefficient  value  of  0.03  was 
used  in  the  computation.  The  results  are  shown  in  Figs.  4.3  to  4.5.  In  these 
figures ,  the  stars  represent  the  case  with  friction  and  the  circles  the  case 
without  friction. 

Figs.  4.3(a)  to  4.3(c)  present  the  normal  forces  acting  on  the  meshing 
gear  teeth  for  the  three  wormgear  positions,  viz.  6*  =  0,  1  and  2  deg.  These 
figures  correspond  to  the  cases  of  Figs.  3.8(a)  to  3.8(c).  The  variation  of 
the  normal  force  on  one  gear  tooth  as  it  moves  through  the  meshing  zone  is 
shown  in  Fig.  4.4,  which  is  essentially  the  three  graphs  of  Figs.  4.3 
combined.  It  is  seen  that  the  difference  between  including  and  not  including 
friction  cannot  be  detected.  This  implies  that,  for  a  given  gear  output 
power,  the  lubricant  film  at  the  worm  and  gear  contact  must  withstand  almost 
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the  same  normal  force  whether  friction  is  present  or  not. 


The  tangential  forces  acting  on  the  worm  thread  corresponding  to  the 
three  wormgear  positions  of  =  0,  1  and  2  deg  are  shown  in  Figs.  4.5(a)  to 
4.5(c).  Here,  including  or  not  including  friction  makes  a  big  difference. 
This  is  because  for  the  selected  wormgear  the  direction  of  the  relative 
velocity,  which  is  the  direction  of  friction,  is  very  close  to  that  of  the 
tangential  force  on  the  worm.  Thus,  a  large  portion  of  the  input  torque  must 
be  spent  to  overcome  the  frictional  torque.  Hence,  effective  lubrication  at 
the  contact  can  drastically  reduce  the  input  power  required  of  the  worm. 


5  Description  of  Computer  Programs 

Two  computer  programs  were  developed  in  this  work,  the  graphics  program 
and  the  analyses  program.  The  graphics  program  is  for  creating  the  3-D 
graphics  of  the  selected  wormgear.  The  analyses  program  is  for  implementing 
the  geometric ,  kinematic  and  load  sharing  analyses . 

5 . 1  Graphics  Program 

The  graphics  program  is  named  NASA.PRG  and  is  written  in  I -DBAS 
programming  language,  IDEAL. 

5.1.1  I-DEAS  Software 

I -DBAS  (Integrated  Design  Engineering  Analysis  System)  is  a  widely  used 
commercial  software  (SDRC,  1991).  It  has  three  program  families:  solid 
modeling,  finite  element  and  advanced  technique.  Solid  modeling  is  used  to 
geometrically  create  an  object.  The  tasks  in  the  solid  modeling  family 
include  construction  geometry,  object  modeling,  mechanism  design,  drawing 
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layout ,  etc . 


Construction  geometry  is  used  to  create  profiles,  wireframes  and  skin 
groups.  A  profile  is  defined  by  a  boundary  that  lies  in  a  plane.  The 
boundary  can  consist  of  lines,  arcs,  circles  and  splines,  which  can  be  defined 
by  points,  centers  and  radii.  In  the  boundary  holes  can  be  defined. 

Object  modeling  is  used  to  design  mechanical  parts  which  are  modeled  as 
objects.  There  are  several  ways  to  create  an  object:  (1)  Use  primitives  such 
as  blocks,  cylinders,  cones  and  spheres.  (2)  Extrude  or  revolve  a  profile 
into  a  solid.  (3)  Create  a  group  of  profiles  and  then  generate  a  solid 
through  the  group.  This  task  also  provides  jointing,  cutting  and  assembling 
options.  Once  two  objects  are  created,  a  third  object  can  be  created  from 
them  by  using  these  options.  The  created  object  can  be  displayed  with  options 
such  as  scale,  hidden  line,  shade  display,  etc.  The  created  object  can  also 
be  rotated,  moved  and  oriented  to  the  required  position;  and  saved  and 
printed. 

I-DEAS  can  be  run  interactively  by  selecting  from  the  menu  on  the  screen 
using  the  mouse,  or  by  inputing  from  the  keyboard  or  through  a  program. 

I-DEAS  has  a  high  level  programming  capability  using  the  language  IDEAL. 
The  commands  can  be  written  in  a  program  file  which  is  then  executed.  The 
programming  capability  permits  arithmetic  expressions,  variables,  mathematic 
functions,  and  algebraic  operations;  it  also  provides  commands  for  branching 
and  do -looping,  as  well  as  controlling  the  input  and  output. 

5.1.2  Description  of  the  Graphics  Program 

The  process  of  creating  a  double -enveloping  wormgear  by  using  I-DEAS  is 
similar  to  that  of  manufacturing  the  wormgear.  This  is  carried  out  as 
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follows : 


First  a  worm  blank  and  a  disc  grinding  wheel  are  created.  Once  these 
two  objects  are  created,  the  worm  can  be  generated  by  cutting  the  worm  blank 
with  the  grinding  wheel .  In  the  worm  creating  process  the  motion  of  the  worm 
blank  and  the  grinding  wheel  is  similar  to  that  of  the  first  enveloping 
process . 

Then  the  gear  blank  is  created.  The  gear  is  created  by  cutting  the  gear 
blank  with  the  previously  created  worm.  The  motion  of  the  gear  blank  and  the 
worm  is  similar  to  that  of  the  second  enveloping  process. 

The  motions  are  simulated  by  rotating  the  blanks  and  cutters  about  their 
respective  axes.  The  rotation  is  divided  into  incremental  steps.  At  each 
step  the  cutting  operation  takes  place .  The  incremental  step  determines  how 
smooth  the  generated  worm  and  gear  surfaces  are.  As  a  result  there  is  a 
difference  between  the  I -DBAS  generated  surface  and  the  manufactured  surface. 
Because  of  the  computer  storage  limitation,  the  steps  can  not  be  chosen  too 
small . 


The  structure  of  the  graphics  program  is  shown  by  the  flow  chart  in  Fig. 

5.1,  which  consists  of  the  following  steps: 

1.  Input  wormgear  dimensions  and  parameters  which  are  needed  in  its 
manufacture . 

2.  Input  the  incremental  stepsizes  ('dtw'  for  generating  the  worm  and  'dtg' 
for  generating  the  gear)  and  the  total  number  of  steps  in  the  cutting 
processes  ('wn'  for  generating  the  worm  and  'gn'  for  generating  the 
gear) . 

3.  Construct  the  grinding  wheels  for  generating  the  worm.  Four  are 
constructed  and  jointed  in  sequence  to  reduce  the  CPU  time,  i.e.  four 
coils  of  the  worm  thread  are  created  for  one  revolution  of  the  worm 
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blank  (dtw*wn=360  deg) . 

4.  Create  the  worm  and  gear  blanks. 

5.  Create  the  worm.  When  the  process  is  completed,  the  created  worm  is 
saved  under  the  name  'worm' . 

6.  Create  the  gear.  When  the  process  is  completed,  the  created  gear  is 
saved  under  the  name  'gear'. 

7.  Display  the  created  worm  and  gear. 

The  program  is  listed  in  Appendix  1.  The  program  works  interactively 
and  ends  at  the  object  modeling  screen.  This  allows  the  user  to  see  the 
created  objects  again  before  exiting  I-DEAS. 


The  procedure  to  execute  the  program  is  as  follows:  (1)  Enter  I-DEAS  and 
come  to  the  construction  geometry  task  screen.  (2)  Run  program  in  the 
construction  geometry  screen.  More  detailed  instructions  are  given  in 
Appendix  1. 

Because  the  worm  thread  of  the  double -enveloping  wormgear  is  asymmetric 
about  its  center,  the  number  of  steps  needed  to  create  the  whole  worm  is 
larger  than  that  for  creating  other  machine  parts,  such  as  the  spur  gear  or 
the  cylindrical  worm.  More  coils  of  the  worm  thread  requires  more  steps. 

Every  step  increases  the  number  of  surfaces  on  the  worm  blank.  The  larger  the 
number  of  surfaces,  the  more  CPU  time  and  computer  disc  space  are  needed  to 
calculate  and  save  the  data.  Once  the  capacity  of  the  computer  or  the  I-DEAS 
space  is  full,  the  screen  stops,  or  the  program  crashes.  This  is  why  there 
are  only  four  coils  of  worm  thread  shown  in  Fig.  2.7.  The  program  was 
executed  on  an  IBM  5080  terminal  at  SUNY- Binghamton. 

From  the  above  description  it  can  be  seen  that  (1)  the  computer 
simulation  is  helpful  to  visualize  the  selected  wormgear,  and  (2)  the  computer 
simulation  can  not  replace  the  analysis,  because  the  numerical  data  of  the 
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generated  surfaces  are  not  accurate  enough  for  the  projected  lubrication 
analysis . 

5.2  Analyses  Program 

5.2.1  Overall  view  of  the  Program 

The  analyses  program  is  named  NASA. FOR  and  is  written  in  FORTRAN.  This 
program  consists  of  a  main  program  and  several  subroutines.  The  main  program 
is  the  driver  program  of  all  the  subroutines  and  may  be  tailored  and  adapted 
to  the  final  design  package.  The  subroutines  fall  in  two  categories:  those 
driven  by  the  main  program  (DIMEN,  WORM,  WORML,  GEAR,  CONTACT,  PITCH  and  LOAD) 
and  those  driven  by  other  subroutines  (EQW,  EQWL,  EQG,  FIRST,  SECOND,  WORMT , 
ROOT  and  NORMAL).  Most  of  the  logical  judgements  are  done  in  the  subroutines 
called  by  the  main  program.  The  evaluations  of  analytical  expressions  are 
done  in  the  subroutines  called  by  other  subroutines . 

The  structure  of  the  analyses  program  is  shown  by  the  flow  chart  in  Fig. 
5.2,  which  consists  of  the  following  steps: 

1.  Calculate  the  wormgear  dimensions  (DIMEN). 

2.  Calculate  the  coordinates  in  S^  of  the  worm  surface  and  the  curvatures 
of  the  worm  surface  (WORM  and  WORML) . 

3.  Calculate  the  coordinates  in  83  of  the  points  on  the  gear  surface  along 
the  contact  line  and  the  curvatures  of  the  gear  surface  (GEAR) . 

4.  Calculate  the  contact  lines  in  S„  on  every  pair  of  the  meshing  teeth  for 
a  given  position  (CONTACT) . 

5.  Calculate  the  coordinates  in  Sj  and  the  kinematic  properties  of  the  worm 
center  points  for  a  given  position  (PITCH) . 

6.  Calculate  the  load  distribution  on  the  meshing  teeth  for  a  given 
position  (LOAD) . 
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Input  of  this  program,  under  the  'parameter'  command,  are  at  the 
beginning  of  the  main  program  and  are  explained  there.  In  the  main  program 
the  given  position  6*  is  named  'tn2' ,  which  should  be  the  input  for  the 
subroutines  CONTACT,  PITCH  and  LOAD. 

Input  to  a  subroutine  are  explained  at  the  beginning  of  the  subroutine. 
The  output  to  be  used  by  other  subroutines  are  described  in  the  subroutine 
statements.  The  output  not  to  be  used  by  other  subroutines  are  stored  in  data 
files  (DIMEN,  WORM,  WORML,  GEAR,  CONTC,  PITCH  and  LOAD). 

Error  flags  would  be  shown  on  the  screen  such  as  'dimension  is  not  big 
enough',  'the  stepsize  is  too  large',  'the  worm  surface  intersects  the  edge 
curve' ,  etc. 

5.2.2  Detailed  Description  of  the  Program 

Main  Program 

The  input  under  the  'parameters'  command  include  the  center  distances  ag 
and  a,  the  speed  ratios  ig  and  i,  the  number  of  worm  thread  Nj,  the  number  of 
meshing  teeth  Np,  the  worm  pitch  diameter  d^,  the  apex  angle  the  rotation 
speed  of  the  worm,  the  output  power  of  the  wormgear,  and  the  friction 
coefficient  at  the  contact. 

The  stepsizes  and  the  dimension  arrays  are  input  under  the  'data' 
command.  The  input  stepsizes  are:  'df'  (for  the  running  parameter  Vi) .  'dt' 
(for  the  running  parameter  62),  'dx'  (in  the  x  axis),  and  'dz'  (in  the  z 
axis).  There  are  four  one -dimensional  common  arrays  with  dimension  nn=600. 


The  given  position  B*  is  named  'tn2'.  The  wormgear  position  of  fi*=0  is 
shown  in  Fig.  3.6.  'tn2'  is  not  to  be  larger  than  the  angle  between  two 
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neighboring  gear  teeth.  The  two  angles  't20'  and  't21'  specify  the  interval 
of  the  running  parameter  62,  which  is  chosen  to  be  slightly  larger  than  the 
meshing  zone. 

The  convergence  criteria  for  root  search  are  specified  by  'ep'  (for 
angles)  and  'epp'  (for  distances).  A  logical  judgement  for  the  number  of 
meshing  teeth  is  set  up  in  the  main  program.  This  is  because  the  actual 
number  of  meshing  teeth  obtained  from  subroutines  CONTACT  and  PITCH  should 
agree,  otherwise  the  results  can  not  be  used  in  the  subroutine  LOAD.  A  flow 
chart  for  the  main  program  is  shown  in  Fig.  5.3. 

Subroutine  DIMEN 


This  subroutine  must  be  run  before  any  other  subroutine,  because  it 
provides  the  wormgear  parameters  and  dimensions.  The  formulas  used  for  the 
calculation  can  be  found  in  AGMA  (1965),  Dudley  (1984),  and  Shigley  and 
Mitchell  (1983). 

The  meshing  zone  is  specified  by  'start'  and  'end',  which  are  the  angles 
where  a  gear  tooth  runs  into  and  out  of  the  meshing  zone  respectively.  'f20' 
and  'f21'  are  the  two  angles  specifying  the  interval  of  the  running  parameter 
^2-  'tao'  is  the  angle  between  two  neighboring  gear  teeth.  The  half  gear 
face  angle,  named  'dlt',  is  used  to  specify  the  interval  in  the  calculation 
of  the  gear  surface.  'reps'  is  an  angle  for  determining  the  position  of  the 
original  generating  plane  that  grinds  the  other  side  of  the  worm  thread.  Fig. 
5.4. 

Subroutine  WORM  and  WORML 


These  subroutines  implement  the  analysis  described  in  Sections  3.3  and 
3.5.  The  two  sides  of  the  worm  thread  are  calculated  by  WORM  and  WORML, 
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respectively . 


Three  logical  judgements  are  set  up.  The  first  checks  if  the  worm 
surface  intersects  the  edge  curve.  The  second  checks  which  branch  of  the 
envelope  surface  is  the  real  worm  surface.  Fig.  3.2.  The  third  examines  the 
tooth  domain. 

Subroutines  EQW  and  EQWL  are  called  by  subroutines  WORM  and  WORML, 
respectively.  A  flow  chart  for  WORM  (and  WORML)  is  shown  in  Fig.  5.5. 

Subroutine  GEAR 


This  subroutine  implements  the  analysis  described  in  Sections  3.4  and 
3.5.  Calculations  for  both  the  conventional  form  and  the  modified  form  are 
provided. 

A  0  value,  named  'tl',  is  first  given.  Then  V*.  named  'fl',  is  varied 
near  the  6  value  to  calculate  the  coordinates  of  the  points  on  the  gear 
surface.  This  is  because  if  V  is  far  away  from  6,  it  is  also  off  the  tooth 
domain.  The  range  of  variation  is  twice  the  gear  face  angle.  This  range  is 
much  less  than  the  range  of  the  running  parameter  ^  in  the  subroutine  WORM. 
Hence,  the  stepsize  *df*  ,  used  here  is  a  tenth  of  the  stepsize  given  in 

the  main  program. 

The  logical  judgements  set  up  in  this  subroutine  consist  of  the  tooth 
domain  examinations  and  sorting  out  the  fictitious  points.  A  flow  chart  for 
GEAR  is  shown  in  Fig.  5.6. 

Subroutine  CONTACT 

This  subroutine  implements  the  analysis  of  Section  3.6.1.  Calculations 
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for  both  the  conventional  form  and  the  modified  form  are  provided.  At  a  given 
position  'tn2',  the  contact  lines  on  all  the  meshing  teeth  are  calculated. 

The  angles  'start'  and  'end'  are  used  in  this  subroutine  to  find  out  how  many 
pairs  of  meshing  teeth  actually  exist  in  the  meshing  zone. 

Similar  to  the  case  of  GEAR,  ip  is  varied  near  the  6  value.  To  obtain 
accurately  the  length  of  the  contact  lines ,  as  many  points  are  needed  as  is 
practicable.  Hence  the  stepsizes  of  ip  and  z  are  chosen  to  be  df/100  and 
dz/10,  where  'df'  and  'dz'  are  the  input  stepsizes  from  the  main  program. 

Four  temporary  one - dimens ional  arrays  are  used  to  save  the  coordinates 
of  the  points  on  the  same  contact  line.  'nn'  gives  the  dimension  of  these 
arrays,  which  is  controlled  by  the  main  program.  If  the  saved  points  exceed 
'nn' ,  an  error  flag  would  appear  on  the  computer  screen.  The  number  'nc' 
(output  name  'ncc')  records  how  many  pairs  of  meshing  teeth  exist  in  the 
meshing  zone  for  the  given  position,  which  should  be  the  same  as  that  obtained 
from  the  subroutine  PITCH. 

Output  of  the  contact  line  length  is  stored  in  the  array  'CONTC',  which 
has  dimension  Np  and  is  the  input  of  the  subroutine  LOAD.  A  flow  chart  for 
CONTACT  is  shown  in  Fig.  5.7. 

Subroutine  PITCH 


This  subroutine  implements  the  analysis  of  Sections  3.4  and  3.6.2.  At 
the  given  position  'tn2'  (named  'fn2'  here),  the  worm  center  point  on  each 
pair  of  the  meshing  teeth  are  obtained  by  first  finding  the  points  on  the  worm 
profile  in  the  mid-plane,  i.e.  y^  =  0,  then  sorting  out  the  center  point  on 
the  profile  by  using  the  subroutine  ROOT.  The  angles  'start'  and  'end'  are 
used  to  check  if  a  worm  coil  falls  in  the  meshing  zone . 
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To  obtain  accurately  the  center  point  coordinates  while  still  saving 
the  CPU  time,  the  stepsize  is  chosen  to  be  df/100,  where  'df'  is  the  input 
stepsize  from  the  main  program.  The  range  of  variation  is  90  deg  around  the 
given  position  'fn2'.  The  index  mw=l  would  mean  that  the  worm  profile 
intersects  the  edge  curve.  The  coordinates  of  the  center  points  obtained  are 
in  Sj.  These  coordinates  and  the  normal  directions  at  the  center  points  are 
stored  in  one -dimensional  arrays  with  the  maximum  dimension  Np,  and  are  input 
for  the  subroutine  LOAD. 

The  ntimber  'nc'  (output  name  'ncp')  records  how  many  pairs  of  meshing 
teeth  exist  in  the  meshing  zone  for  the  given  position,  which  should  be  the 
same  as  that  obtained  from  the  subroutine  CONTACT.  A  flow  chart  for  PITCH  is 
shown  in  Fig .  5.8. 

Subroutine  LOAD 

This  subroutine  implements  the  analysis  of  Sections  4. 
subroutine  is  run,  the  subroutines  PITCH  and  CONTACT  must  be 
provide  the  input  for  this  subroutine. 

Subroutine  ROOT 

This  subroutine  searchs  the  root  in  a  given  interval  where  only  one  root 
exists.  The  root  searching  is  based  on  the  idea  of  bisection  (Press,  1986), 
but  three  sections  are  used  instead  to  reduce  the  CPU  time.  Here  the  root  is 
a  point  at  which  the  radius  of  the  worm  is  equal  to  the  radius  of  the  worm 
pitch  cylinder  at  the  same  axial  position.  The  criteria  used  to  search  the 
center  point  are  the  difference  of  the  two  radii  or  the  angles  between  the  two 
grid  points. 


Before  this 
run  so  as  to 
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6  Summary 


Based  on  the  consideration  of  achieving  better  contact  properties 
(longer  contact  lines  and  a  larger  angle  between  the  contact  line  and  the 
relative  sliding  velocity),  a  double -enveloping  wormgear  was  selected.  The 
worm  was  generated  by  a  disk  grinding  wheel.  The  center  distance  and  the 
speed  ratio  were  respectively  the  same  in  the  first  and  the  second  enveloping 
processes  (conventional  form) . 

Based  on  the  suggested  speed  ratio  and  center  distance  for  helicopter 
applications,  the  other  dimensions  and  geometric  parameters  of  the  wormgear 
were  laid  out  basically  in  accordance  with  the  AGMA  gear  standards.  However, 
some  of  these  quantities  were  modified  as  a  result  of  compromising  between 
conflicting  requirements.  It  is  felt  that  for  a  novel  application  like  ours 
that  is  not  covered  by  past  experiences,  some  deviations  from  the  common 
practice  would  be  inevitable. 

3-D  graphics  representation  of  the  selected  wormgear  was  developed  using 
the  I-DEAS  software.  The  software  works  like  a  cutting  tool  generating  a 
surface  from  a  blank,  but  in  discrete  steps.  Because  of  the  computer  storage 
limitation,  the  steps  could  not  be  made  too  small.  As  a  result,  the  numerical 
data  of  the  generated  surfaces  were  not  accurate  enough  for  the  projected 
lubrication  analysis.  However,  the  graphics  representation  does  convey  a 
clear  notion  of  the  proportion  of  the  selected  wormgear. 

The  kinematics  of  meshing  was  analyzed  with  the  method  of  enveloping. 

The  method  is  simple  and  elegant,  and  can  be  used  to  construct  conjugate 
surfaces  in  other  contacting  mechanisms  as  well.  Using  this  method  the  shape 
of  the  worm  and  gear  surfaces,  the  contact  lines,  the  normal  directions,  the 
relative  normal  curvatures ,  and  the  relative  sliding  velocities  at  the  contact 
points  were  obtained.  With  these  results  the  clearance  function  between  the 
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mesliing  surfaces  can  also  be  obtained,  which  is  needed  for  the  projected 
lubrication  analysis. 

Based  on  the  contact  line  pattern  revealed  by  the  kinematics  analysis,  a 
hypothesis  of  load  sharing  was  made.  Using  this  hypothesis  the  variation  of 
the  normal  contact  force  acting  on  a  gear  tooth  as  it  ran  through  the  meshing 
zone,  as  well  as  the  distribution  of  the  normal  forces  on  all  the  meshing  gear 
teeth  at  any  instant,  were  obtained.  This  information  is  needed  in  the 
projected  lubrication  analysis.  The  required  torque  output  of  the  worm  was 
likewise  obtained.  The  results  clearly  demonstrated  the  advantage  and  need  of 
implementing  effective  lubrication  at  the  contact. 

A  computer  program,  written  in  FORTRAN,  was  developed  to  carry  out  the 
computations.  The  structure  of  the  computer  program  and  its  usage  are 
described  in  detail. 
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Table  1  Dimensions  of  the  Selected  Wormgear 


TERM  SYMBOL 

VALUE 

FORMTITA  BASED 

REFERENCE 

speed  ratio 

i 

110 

given 

normal  pressure 
angle 

“n 

20  deg 

Dudley,  1984 
p.  3.71 

center  distance 

a 

880  mm 

selected 

number  of  worm 
threads 

Ni 

1 

selected 

number  of  gear 
teeth 

N2 

110 

N2  =  Ni*i 

pitch  diameter 
of  worm* 

di 

300  mm 

di  =  a°®”Ad  (E) 

(kd  =  1.7-2. 2) 

Dudley,  1984 
p.  3.70 

pitch  diameter 
of  gear 

^2 

1460  mm 

dg  =  2*a-di 

axial  pitch 

Px 

41.6975  mm 

p,,  =  9r*d2/N2 

Dudley,  1984 
Table  3.32 

worm  lead  angle* 

7 

2.533  deg 

tan7  =  Px*Ni/jr/di 

Dudley,  1984 
Eq.  3.38 

normal  cirlular 
pitch 

Pn 

41.6568  mm 

Pn  =  Px*COS7 

Dudley ,  1984 
Table  3.32 

axial  pressure 
angle* 

20.018  deg 

tana^  =  tanci:jj/cos7 

Dudley ,  1984 
Eq.  3.40 

Aa 

0.818  deg 

=  Px/2/d2 

Dudley,  1984 
p.  3.72 

base  circle 
diameter 

db 

519.3182  mm 

db  =  d2*sin(aA^o:) 

Dudley,  1984 
Eq.  3.41 

module 

Xn. 

13.2727 

mm/tooth 

Xn,  =  PxA 

Dudley,  1984 
Table  3.32 

number  of 
meshing  teeth 

Np 

10 

selected 

half  angle  of 
meshing 

V-a 

15.627  deg 

V>a  =  7r*(Np-0.45)/N2 

Shen,  1983 
Table  10.1 

start  angle 
of  meshing 

V'f 

5.209  deg 

V>f  =  sin'^(dt,/d2)-^a 

Shen,  1983 
Table  10.1 

thickness  of 
worm  thread 

Si 

18.7639  mm 

Si  =  0.45*p^ 

Dudley,  1984 
p.  3.73 

thickness  of 
gear  tooth 

S2 

22.9336  mm 

S2  =  0.55*P3, 

Dudley,  1984 
P.  3.73 
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whole  depth 

ht 

20.8284  mm 

working  depth 

bic 

18.7455  mm 

addendum 

ba 

9.3728  mm 

clearance 

c 

2.0829  mm 

throat  diameter 
of  worm* 

dal 

318.7455  mm 

root  diameter 
of  worm* 

dfi 

277.0888  mm 

throat  diameter 
of  gear 

da2 

1478.7455  mm 

root  diameter 
of  gear 

df2 

1437.0888  mm 

face  width  of 
worm 

bi 

393  mm 

face  width  of 
gear 

^2 

249  mm 

apex  angle  of 
generating  plane 

P 

5  deg 

bfc  =  0.5*p„ 

Dudley,  1984 
Table  3.31 

bk  =  0.45*p„ 

Dudley,  1984 
Table  3.31 

b,  =  0.225*p, 

c  =  bt-bj. 

Dudley,  1984 
Table  3.31 

dal  =  di+2*b. 

Shen,  1983 
Table  10.1 

dfi  =  dai-2*bt 

Dudley,  1984 
Eq.  3.39 

da2  ~  d2+2*b3 

Dudley,  1984 
Table  3.32 

df2  “  da2-2*bt 

Shen,  1983 
Table  10.1 

bi  =  d2*sinV'a 

Shen,  1983 
Table  10.1 

b2  =  (0.9-1.0)*dfi 

Shen,  1983 
Table  10.1 

selected 

*  At  the  center  o£  worm. 
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Figures 


Fig.  2.1 
Fig.  2.2 
Fig.  2.3 
Fig.  2.4 
Fig.  2.5 
Fig.  2.6 
Fig.  2.7 
Fig.  2.8 
Fig.  2.9 
Fig.  3.1 
Fig.  3.2 
Fig.  3.3 
Fig.  3.4 
Fig.  3.5 
Fig.  3.6 
Fig.  3.7 
Fig.  3.8 

Fig.  3.9 
Fig.  3.10 
Fig.  4.1 
Fig.  4.2 
Fig.  4.3 

Fig.  4.4 
Fig.  4.5 

Fig.  5.1 


Single -enveloping  wormgear  (cylindrical  wormgear) 

Double -enveloping  wormgear 

Worm  ground  by  a  disc  grinding  wheel 

Worm  ground  by  a  cone  grinding  wheel 

Double -enveloping  wormgear  layout  and  nomenclature 

Selected  wormgear 

Selected  worm 

Selected  gear 

Close-up  view  of  the  selected  gear  teeth 
Coordinate  systems 

Envelope  surface  (worm)  and  its  edge  curve 

Contact  lines  on  a  gear  tooth  of  the  modified  form 

Contact  lines  on  a  gear  tooth  of  the  conventional  form 

Worm  and  gear  tooth  domains 

Wormgear  position  at  5*  =  0 

Contact  lines  on  the  meshing  gear  teeth 

Length  of  contact  lines  on  the  meshing  teeth 

(a)  position  of  d*j=0  deg 

(b)  position  of  ^*=1  deg 

(c)  position  of  6  =2  deg 

Variation  of  contact  line  length  on  a  gear  tooth 
Radial  position  of  the  edge  curve 
Normal  force  and  its  components 
Friction  force  and  its  components 

Normal  forces  on  the  gear  teeth 

(a)  position  of  0*=O  deg 

(b)  position  of  0  =1  deg 

(c)  position  of  d*=2  deg 

Variation  of  normal  force  on  a  gear  tooth 

Tangential  forces  on  the  worm 

(a)  position  of  fi*=0  deg 

(b)  position  of  0*=1  deg 

(c)  position  of  6  =2  deg 

Structure  of  graphics  program 
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Fig.  5.2 
Fig.  5.3 

Fig.  5.4 

Fig.  5.5 
Fig.  5.6 
Fig.  5.7 
Fig.  5.8 


Structure  of  analyses  program 
Flow  chart  of  main  program 

Two  sides  of  the  worm  thread 


Flow  chart 
Flow  chart 
Flow  chart 
Flow  chart 


of  subroutine  UORM  (also  WORML) 
of  subroutine  GEAR 
of  subroutine  CONTACT 
of  subroutine  PITCH 
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Fig.  2.2  Double-enveloping  wormgear 


46 


Fig,  2.4  Worm  ground  by  a  cone  grinding 
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Fig.  2.5  Double-enveloping  wormgear  layout  and  nomenclature 


Fig.  2.6  Selected  wormgear 
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Fig.  2.8  Selected  gear 


Y? 
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Fig.  2.9  Close-up  view  of  the  selected  gear  teeth 
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Fig.  3.1  Coordinate  systems 
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Fig.  3.2  Envelope  surface  (worm)  and  its  edge  curve 
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Fig.  3.7  Contact  lines  on  the  meshing  gear  teeth  (b) 


125.00 


No.8 


o  o  o  o  o  gear 

*  *  +  *  *  worm 


boundary 


)0  ®  OB’*  *  *  +  * 


JOKDCxa  +  * 


>  ®=OO’I0+  * 


r-the  second  enveloping  contact  line 


-the  first  enveloping  contact  line 


>€>€>=G* 


>09  9  * 


>0€)l'Odl?  *  *  * 


125.00  I  I  M  I  I  I  1  M  I  I  i  I  I  I  1  1  i  t  I  i  1  I  I  I  i  i  M  I 

600.00  650.00  700.00  750.00  800.00 


Fig.  3.7  Contact  lines  on  the  meshing  gear  teeth  (d) 
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Fig.  3.7  Contact  lines  on  the  meshing  gear  teeth  (f) 
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Fig.  3.8  Length  of  contact  lines  on  the  meshing  teeth 
(a)  position  of  e*=0  deg 
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Fig.  3.8  Length  of  contact  lines  on  the  meshing  teeth 
(c)  position  of  0*=2  deg 
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Fig.  3.10  Radial  position  of  the  edge  curve 


Fig.  4.1  Normal  force  and  its  components 
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Fig,  4.2  Friction  force  and  its  components 


72 


Normal  Force  (N) 


70000.00 


60000.00  H 


50000.00  H 


40000.00  H 


30000.00  H 


position  0*  =  0  (deg) 


with  friction 
o  &€)  GO  without  friction 


©-* 


\ 

\ 

\ 


\ 


0 


20000.00  — I  I  I  I  I  I  I  I  I  I  I  I  I 'I  I  I  I  I  I  IT- 


Ill  I  I  1 1 1'  I  I  I  I  'T"!  i  1  I  I  I  I  r 

0.00  10.00  20-00  30.00  40.00 


Tooth  Position  0t  (deg) 


Fig.  4.3  Normal  forces  on  the  gear  teeth 
(a)  position  of  e*=0  deg 
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Fig.  4.3  Normal  forces  on  the  gear  teeth 
(b)  position  of  e*=1  deg 
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Fig.  4.4  Variation  of  normal  force  on  a  gear  tooth 
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Fig.  4.5  Tangential  forces^  on  the  worm 
(a)  position  of  0*=O  deg 
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Fig.  4.5  Tangential  forces  on  the  worm 
(b)  position  of  6*=1  deg 
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Fig.  4.5  Tangential  forces  on  the  worm 
(c)  position  of  0*=2  deg 


Fig.  5.1  Structure  of  graphics  program 
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MAIN 


Fig.  5.3  Flow  chart  of  main  program 
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Fig.  5.5  Flow  chart  of  subroutine  WORM  (also  WORML) 


CaU  SECOND 


Fig.  5.7  Flow  chart  of  subroutine  CONTACT 
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Rg.  5.7  (continued) 


Fig.  5.8  Flow  chart  of  subroutine  PITCH 
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CaU  NORMAL 
1  calculate  normal  direction 
and  relative  velodty. 


Fig.  5.8  (continued) 
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Appendix  1  Graphics  Program  (NASA.PRG) 


The  procedure  to  execute  'NASA.PRG* 

(1)  From  the  SUNY- Binghamton  operating  system  prompt  (assumed  here  to  be  $) 
$  ready;  type  'runideas'  <retum> 

$  jtunp  to  IBM  screen  (use  Alt-jump) 

$  user  i.d.:  type  'ideas2' 

$  password:  type  'no' 

$  jump  back  to  BINGVMB  screen  tjrpe  'none'  (use  Alt -jump) 

$  disk  space;  type  a  number  such  as  '15000' 

(waiting  for  temporary  disk  to  be  created) 

$  user  i.d.:  t3rpe  your  own  i.d.  ntimber  and  read  password 
$  terminal  type:  type  'i5080' 

$  jump  screen  (use  Alt-jump) 

(waiting  for  I -DBAS  screen  to  appear) 

(2)  Under  the  I-DEAS  screen,  from  the  I-DEAS  operating  prompt  (assumed  here 
to  be  #) : 

Enter  model  file  name  (none)//  <retum> 

New  model  file?  (yes)#  <return> 

Enter  model  description  (none)#  <return> 

Enter  system  of  units  (SI)#  <return> 

Enter  operation  mode  (Normal_Operation)#  <retum> 

(3)  Then  choose  which  "family"  of  I-DEAS  you  want.  Notice  that  choices  are 
in  the  form  of  a  'menu' .  So  you  can  pick  a  choice  with  the  mouse  or  type 
under  the  I-DEAS  operating  prompt  (assumed  here  to  be  #) : 

Select  menu#  SO  <return> 

(waiting  for  Solid  Modeling  Screen  to  appear) 

Select  menu//  CG  <return> 
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Select  menu#  MF  <return> 

Select  menu#  PR  <return> 

Select  menu#  R  <return> 

File  name#  NASA  <return> 

Then  the  program  is  run  interactively.  All  the  results  will  be  seen  from  the 
screen. 

(4)  The  created  object  can  be  observed  again  before  exiting  I-DEAS  system. 
The  procedure  to  observe  the  created  object  is: 

Select  menu#  MA  <return> 

Select  menu#  G  <return> 

Enter  object  name#  type  object  name  (worm  or  gear) 

(the  created  object  appears  on  the  screen) 

Once  the  created  object  is  obtained,  choose  from  the  menu  options  such  as 
orient,  view,  display  and  AU  (auto  scale)  to  observe  the  object. 

(5)  When  the  observation  is  completed,  put  the  created  object  away  from  the 
screen  by: 

Select  menu#  /  <return> 

Select  menu#  MA  <return> 

Select  menu#  PU  <return> 

(6)  When  all  the  observations  are  completed,  if  needed  the  created  object 
can  be  saved  by: 

Select  menu#  /  <return> 

Select  menu#  MA  <return> 

Select  menu#  WR  <return> 

Select  menu#  U  <retum> 

Enter  the  file  name#  type  the  name  you  choose  <return> 
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(7)  Exit  I-DEAS  screen  by: 

Select  menu#  /  <retum> 
Select  menu#  type  'exit' 


<return> 
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Program  Llstlnt 


****************************************************************** 
This  is  an  ideas  program  to  construct  a  wormgear  set  by  the 
principle  of  enevlop.  The  constructing  process  is  same  as 
the  manufacturing  process  exactly.  To  cut  four  worm  threads 
the  'wn'  must  be  equal  to  360/dtw  and  'gn'  must  be  4*360/dtg. 


The  basic  parameters  provided  are  : 


aO . center  distance  of  cutting  (a( 

db . base  circle  diameter  (dl 

bl - face  width  of  worm  (b, 

hi - face  width  of  gear  (b: 

d2 - pitch  diameter  of  gear  (d: 

da2 - thread  diameter  of  gear  (d. 

df2----root  diameter  of  gear  (d 

dfl - root  diameter  of  worm  (d 

ratio --speed  ratio  of  wormgear  (r 

n2 - nimiber  of  gear  teeth  (n 

dtl--- -angle  of  worm  turning  in  each  cut 
dt2--- -angle  of  gear  turning  in  each  cut 

wn . times  of  worm  cutting 

gn . times  of  gear  cutting  _ 


(a0=880.0  mm) 
(db-519.3182  mm) 
(bl=383.0  mm) 
(b2-249.0  mm) 
(d2=1460.00  mm) 
(da2=1478.7455  mm) 
(df 2=143 7. 08 8 8  mm) 
(df 1=277. 0888  mm) 
(ratio=110) 
(n2=110) 


gn . times  or  gear  cuccing 

bta - apex  angle  of  generating  plane  (bta=10  deg) 

psif angle  of  meshing  start  (psif=5.209  deg) 

reps angle  of  two  cutting  surfaces  (reps=39 . 8724  deg) 

**  'reps'  is  calculated  by  the  following  formular: 
reps=180-{2*[90-arcsin(db/d2)]+s2/r2}  ,  here 
s2=0.55*pi*d2/n2  i.e.  s2/r2=l . I*180/n2  (deg) 
****************************************************************** 


=  basic  parameters  input 

II  ^11 

#echo  list 

#a0=880.0 

#n2=110 

#ratio=110 

#db=519 . 3182 

#bta=5 . 00 

#dl=300.0 

#d2=1460.0 

#dal=318.7455 

#da2=1478 . 7455 

#dfl=277.0888 

#df2=1437.0888 

#bl=393.0 

#b2=249.0 

#rb=db/2.0 
#rl=dl/2.0 
#r2=d2/2.0 
#ral=dal/2.0 
#ra2=da2/2.0 
#rfl=df 1/2.0 
#rf2=df 2/2.0 
#dh=2.0829 
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K  :  #psif=5.209 
K  :  #reps-39.8724 

K  :  #aa-a0-sqrt((rf2+dh)**2-bl*bl/4.0) 

K  :  #bb=sqrt(ra2*ra2-rb*rb) 

K  :  #cc=a0-sqrt(rl*rl-b2*b2/4.0) 

K  :  #dd=a0-sqrt(ra2*ra2-bl*bl/4.0) 

K  :  #ee=aO-ra2 
K  :  #dtw=8 . 00 
K  :  #dtg=15.0 
K  :  #wn=5 
K  ;  #gn=3 
C  ; 

C  :  . — 

C  :  =  creat  2-D  profile  of  worm  blank  = 

C  :  -- 

C  : 

K  :  "CR" 

K  :  "PR" 

K  :  "K" 

K  :  0,r2 
K  :  a0-cc,r2 
K  *  **F** 

K  :  a0-cc,bl/2.0 
K  * 

K  aa,bl/2.0 
K  :  "F" 

K  :  a0-rf2-db,0 
K  ;  aa,-bl/2.0 
K  :  "F" 

K  :  aO-cc, -bl/2.0 
K  :  "F" 

K  :  a0-cc,-r2 
K  :  "F" 

K  :  0,-r2 
K  :  "F" 

K  :  "C" 

K  ; 

K  :  "MA" 

K  :  "PR" 

K  :  "STO" 

K  :  pi 
K  :  "/" 

K  :  "DE" 

K  ;  "P" 

K  :  "Y" 

K  :  "/" 

C  : 

C  :  .  - ^ -  - 

C  :  =  creat  2-D  profile  of  worm  root  blank  = 

C  :  —  ■ 

C  : 

K  :  "CR" 

K  ;  "PR" 

V  •  ”1^” 

K  :  0,bl/2 
K  :  dd,bl/2 
K  :  "F" 

K  :  ee,0 
K  :  dd, -bl/2.0 

.  tl  p  If 

K  :  0, -bl/2.0 
K  :  "F" 
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K  : 

"C" 

K  : 

■1  yii 

K  : 

"MA" 

K  : 

"PR" 

K  : 

"STO" 

K  : 

p2 

K  : 

If  yif 

K  : 

”DE" 

K  : 

It  pit 

K  : 

iiY” 

K  : 
C  : 

II  yii 

C  : 

C  : 

=  creat 

3-D  gear  blank  = 

C  : 
C  : 

K  : 

"TA" 

K  : 

"0" 

K  : 

11  ^11 

K  : 

"CR" 

K  : 

"DEF" 

K  ; 

II CY" 

K  : 

npil 

K  : 

32 

K  : 

II  Dll 

K  : 

"D" 

K  : 

II  CY" 

K  : 

cc,b2 

K  : 

II  Jn 

K  : 

"OR" 

K  : 
K  : 

"RO" 

K  : 

90,0,0 

K  : 

"MO" 

K  : 

"D" 

K  : 

0,1,0 

K  : 

aO 

K  : 

II 

K  : 

"MA" 

K  : 

npu" 

K  : 
K  ; 
C  : 

g2 

C  : 

C  : 

=  creat 

a  3-D  worm  cutter  = 

C  : 
C  : 

K  : 

II  yii 

K  : 

"CR" 

K  ; 

"SPH" 

K  ; 

"A" 

K  : 

"C" 

K  : 

6 

K  : 

"D" 

K  : 

ra2+2*c3h 

K  : 

II  yil 

K  ; 

"CO" 

K  ; 

lip  II 

K  : 

lip  II 

K  : 

"K" 

K  : 

-rb,0,0 

K  ; 

-rb,a0,0 
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0,0,rb/taTi(bta) 


rb*cos(reps) , ■rb*sin(reps) ,0 

rb/cos(reps) ,0,0 

0,0,rb/tan(bta) 


II 

"OR" 

"MO" 

II D" 

0,1,0 

aO 

"R" 

0 ,  aO ,  0 
0,0,psif 

liy  II 

"MA" 

"STO" 

c2 


=  creat  four  connected  3-D  worm  cutters  = 


II  yii 

"OR" 

"Rll 

"K" 

0 ,  aO ,  0 
0,0,360/n2 

II  ^11 

"MA" 

"G" 

c2 


II  yi! 

"MA" 

"STO" 

c2 

iiY** 

II  ^it 

"OR" 

"Rll 

"K" 

0 ,  aO ,  0 

0,0,2*360/n2 

I!  ^11 

"MA" 

"G" 

c2 
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0,0,-b2/2 

"D" 

0.0,1 

ItptI 

tip  II 
II  pit 

0,0,b2/2 

"D" 

0.0,1 

"N" 


"K” 

0 , aa+100 , 0 
"D" 

0,1,0 

"N" 

It  Jn 

"MA" 

"STO" 

c2 

iiY" 
n  Jn 

"MO" 

"DE" 

HY" 

It  yii 

"MA" 

"PU" 

c2 

IIY" 


=  creat  3-D  worm  blank  = 


II  yil 

"CR" 

"DEF" 

"R" 

"X" 

32 

"F" 

"A" 

11 

"D" 

"D" 

"REV" 

pi 


0,0,90 

II  ^11 

"MA" 
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=  3-D  worm  cutting 


K  :  "STO" 

K  :  wl 
K  :  "G" 

K  :  c2 
C  : 

C  : 

C  : 

C  : 

C  : 

K  : 

K  :  "CO" 

K  :  "CU" 

K  : 

K  : 

K  :  #time=l 
K  :  #start: 

K  :  y/tl=time*dtw 
K  :  #t2=tl/ratio 
K  :  "OR" 

K  :  "R" 

K  : 

K  :  - dtw ,0,0 
K  :  "/" 

K  :  "MA" 

K  ;  "G" 

K  :  c2 
K  :  "/" 

K  :  "OR" 

K  :  "R" 

K  * 

K  :  0,a0,0 
K  :  0.0, t2 
K  :  "/" 

K  :  "CO" 

K  :  "CU" 

K  ; 

K  :  "/" 

K  :  #time=time+l 

K  :  #if(time  It  wn)  then  goto  start 
K  :  "/" 

K  :  "OR" 

K  :  "R" 

K  : 

K  :  (time-l)*dtl,0,0 
K  :  "/" 

K  :  "MA" 

K  ;  "STO" 

K  ;  worm 
K  : 

K  :  "/" 

K  :  "MO" 

K  :  "DE" 

K  :  "Y" 

K  :  "/" 

K  :  "CR" 

K  :  "REV" 

K  :  p2 
K  : 

K  :  "/" 

K  :  "OR" 

K  :  "R" 

K  : 
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0,0,90 

11^  If 

"CO" 

IIJII 


worm 

ify** 


=  3-D  gear  cutting  = 


II  ^1! 

"MA" 

II G" 

g2 

II G" 

worm 

II 

"CO" 

"CU" 

II  yll 

y/time=l 

y/again: 

#tl=time*dtg 

#dt2=dtg/ratio 


0 ,  aO ,  0 
0,0,  dt2 

liy  11 

"MA" 

"G" 

worm 

II  ^11 

"OR" 

"RII 

-tl,0,0 

II  ^11 

"CO" 

"CU" 


II  ^11 

#time=time+l 

Hf(time  It  gn)  then  goto  again 


**K” 

0 ,  aO ,  0 

0,0,-(time-l)*dt2 

iiyii 

"MA" 

"STO" 

gear 
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K  :  "MA" 

I ;  "FU" 


IMF 

MY” 

#echo  all 


=  view  worm  = 


"G" 

worm 

II  ^11 
liyil 
•Igft 

1,1,1,60 

"MG" 

II  ^11 

"MG" 

II Y** 

"DO" 

"HI" 

"E" 

,1  yti 

"DE" 

,1 Y” 


=  view  gear  = 


gear 

nyu 

"V" 

"RM" 

0,0,120 

„ 

"DO" 

"SH" 

"E" 

****  end  of  session  **** 
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Appendix  2  Listing  of  the  Analyses  Program  (NASA. FOR) 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


******************************************************************** 
PROGRAM  MAIN: 

This  is  a  part  of  a  developing  software  package  which  is  used  to 
analyse  the  double-enveloping  wormgear  fluid  film  lubrication. 

This  part  includes: 

(1)  calculation  of  the  wormgear  dimensions; 

(2)  calculation  of  the  coordinates  of  worm  and  gear  surfaces 
in  SI  and  S2  systems; 

(3)  calculation  of  the  contact  lines  on  all  pairs  of  meshing 
teeth  for  a  given  position; 

(4)  calculation  of  the  coordinates  and  relative  velocity  at 
the  worm  pitch  points  on  all  pairs  of  meshing  teeth  for 
a  given  position; 

(5)  calculation  of  load  sharing  among  all  meshing  teeth  for 
a  given  position. 


To 


run  this 
aO 
a2 
bta 
dl 
dfl 
dt2 
dx 
dz 

ep , epp 

nl 

np 

nn 

rfn 

ratio 

ratio2 

power 

rpm 

tn2 

uf 


program  the  input  parameters  are: 
center  distance  of  the  first  enveloping 
center  distance  of  wormgear 
apex  angle  of  the  generating  plane 
pitch  diameter  of  worm 

51  running  step  size  in  the  first  enveloping 

52  running  step  size  in  the  second  enveloping 
interval  along  x  axis 

interval  along  z  axis 

convergence  criterion  for  root  searching 

number  of  worm  thread 

number  of  meshing  teeth 

number  of  dimension 

normal  presure  angle 

speed  ratio  of  the  first  enveloping 

speed  ratio  of  wormgear 

output  power  from  gear  shaft 

input  speed  from  worm  shaft 

fiven  gear  position 
riction  coefficient 


(mm) 

(mm) 

(deg) 

(mm) 

(deg) 

(deg) 

(mm) 

(mm) 


(deg) 


(kw) 

(rpm) 

(deg) 


The  variables  in  the  main  program  are: 

t20  ----  running  angle  of  the  second  enveloping  start  (rad) 
t21  running  angle  of  the  second  enveloping  end  (rad) 


Subroutine  s  included : 


(1)  Called  by  MAIN  Program 


DIMEN 

WORM 

WORML 

GEAR 

CONTACT 

PITCH 


develop  the  dimensions  of  wormgear. 

develop  the  worm  surface  in  SI .  Here  the  worm  is 

ground  by  a  disc  grinding  wheel. 

develop  the  worm  surface  on  the  other  side,  in  SI. 
develop  the  gear  surface  in  S2. 

calculate  the  length  of  contact  lines  on  each  pair 

of  meshing  teeth  at  the  given  position. 

finding  the  pitch  points  on  the  worm  surface,  the 
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LOAD 


c 

c 

c 

c 

c 

C 

C 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


c 

c 

c 

c 

c 


c 


normal  directions  and  relative  velocities  on  each 
pair  of  meshing  teeth  for  a  given  position, 
calculate  the  normal  force  on  each  pair  of  meshing 
teeth  and  the  torque  sharing  among  the  worm  meshing 
teeth. 


(2)  Called  by  Subroutine 

E(jG  ttie  gear  surface  eptiou  in  11 


EQWL 

FIRST 

NORMAL 


tie  worm  surface  equation  in  si 

the  worm  surface  equation  of  the  other  side  in  SI 
calculation  of  the  first  enveloping  surface 
calculation  of  the  properties  of  the  worm  pitch  points 


ROOT  finding  the  pitch  point 

SECOND  calculation  of  the  second  enveloping  surface 

WORMT  calculation  of  the  worm  profile  on  mid-plane  in  Sj 


Outputs  of  the  program  are  written  in  the  data  files  under 
the  names  of  the  subroutines .  Units  used  in  the  output  data  are ; 

rad  for  angle 

mm  for  length 

newton  for  force . 


PROGRAM  MAIN 

implicit  double  precision(a-h,o-z) 

parameter (bta=5 . 0 , a0=880 . 00 , a2=880 . 00 , ratio“110 . , ratio2=110 , 
&  rfn=20 , nl=l , np=10 , dl=300 . 0 , epp=l . Oe -  5 , ep=l . Oe -  3 , 

&  power=4000 , rpm=14280 ,uf=0 . 03 ,nn=600) 

dimension  contc(np) , tp2(np) ,fp2(np) ,pxj (np) ,pzj (np) , 

&  xpn(np) ,ypn(np) ,zpn(np) ,vx(np) ,vy(np) ,vz(np) , 

&  wtor(np) ,wtorf (np) ,forn(np) , fornf (np) , 

&  exaa(nn) ,exbb(nn) ,excc(nn) ,exdd(nn) 

data  df 1/10 . 0/ , dt2/l . 0/ , dx/5 . 0/ , dz/5 . 0/ , tn2/0 . 0/ 

open (l,file=' /dimen  dat') 
open(2,file='/worm  dat') 
open(3 , file='/worml  dat') 
open(4, file='/gear  dat') 
open(5 , file='/contc  dat') 
open(7 , file=' /normal  dat') 
open(8 , file='/pitch  dat') 
open(9 , file='/load  dat') 

pi=4 . d0*datan(l . dO) 
bt=bta*pi/180 . dO 
tb=dtan(bt) 
dt=dt2*pi/180.d0 
df-dfl*pi/180.d0 
rfnr=rfn*pi/180 . dO 


calculation  of  the 


limensions 


call  dimen (aO , a2 , ratio2 ,nl ,np , dl , rfnr ,bt , 
&  bl,b2,cc,ddd,dlt,f20,f21, 

&  rb , rl , r2 , ral , ra2 , rf 1 , rf 2 , rmax2 , 

&  reps , tt2 , cmax.cmin, start, end) 


print*,'  DIMEN  done' 
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c 

c 


call  wonn(aO,a2, ratio, rb.bt.tb, 

&  f20,f21,bl,ra2,rf2,cc,df ,dz) 


c 

c 

c 


call  worml (a0,a2, ratio, rb,bt,tb, 

&  f 20 , f21 ,bl , ra2 , rf 2 , cc , df , dz , reps) 

Print*, 'WORM  and  WORML  done' 


calcukation  of  the  gear  surfaces  along  contact  lines  | 


t20=f20 

t21=f21 

call  gear (aO , a2 , ratio , ratio2 , rb , tb , bt , rmax2 , cmax , cmin , 

6e  ra2 ,  rf  2 ,  r2 ,  ral ,  rfl ,  rl ,  cc ,  ddd,bl  ,b2 , 

&  f20,f21,df,t20,t21,dt,dlt,dz) 

print* , ' GEAR  done ' 

I  calculation  of  the  contact  line  at  the  given  position  tn2  [ 


call  contact (aO , a2 , ratio , ratio2 ,np , rb , tb ,bt , dh , ddd, 

&  rl , ral , rfl , r2 , ra2 , rf 2 , rmax2 , cmax , cmin, bl ,b2 , 

&  f20 , f21 , df , t20 , t21, tt2 , tn2 , start , end, dlt , dz , 

&  contc , ncc , nn , exaa , exbb , excc , exdd) 

print*, 'CONTACT  done' 


calculate  the  characters  of  pitch  points  on  the  worm  surrace 
at  the  given  position  'tn2' 


call  pitch (aO , a2 , ratio , ratio2 ,np , rb , tb ,bt , cc , ddd, 

&  rl , ral , rfl , r2 , ra2 , rf 2 , rmax2 , cmax , cmin , bl ,b2 , 

&  f 20 , f 21 , df , t20 , t21 , tt2 , tn2 , start , end , dlt , dz , 

&  nn , exaa , exbb , excc , exdd , ep , epp , 

&  tp2,fp2,pxj ,pzj ,xpn,ypn,zpn,vx,vy,vz,ncp) 

print*, 'PITCH  done' 


I  calculate  the  load  distribution  at  the  given  position  tn2 


if (ncc.ne.ncp)  then 

print*,  '  ncc.ne.ncp  check  parameters' 
print*, ncc, ncp 
endif 

call  load  (a2,ratio2,r2,np,power,rpm,uf ,ncc,tp2,pxj , 

&  pzj  , xpn,ypn,zpn, contc, wtor,wtorf,forn,fomf) 

print*, 'LOAD  done' 

stop 

end 

*****************  end  of  MAIN  PROGRAM  ****************************** 


c 

c 

c 
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c 


Subroutine  DIMEN  called  by  MAIN 

This  subroutine  is  used  to  develop  the  dimensions  of  a  wormgear 


lo  1  this  stoutine  tie  iopt  piieters  are; 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


a6  center  distance  of  tne  first  enveloping 

a2  center  distance  of  wormgear 

bta  apex  angle  of  the  generating  plane 

dl  pitch  diameter  of  worm 

nl  number  of  worm  thread 

nbao  number  of  meshing  teeth 

ratio  speed  ratio  of  wormgear 

rfn  normal  presure  angle 


(mm) 

(mm) 

(deg) 

(mm) 


(deg) 


The  output  of  the  subroutine  are: 

bl  face  width  of  worm  (mm) 

b2  face  width  of  gear  (mm) 

cmax  maximum  thread  radius  of  the  gear  hob  (mm) 

cmin  minimum  root  radius  of  the  gear  hob  (mm) 

dh  clearance  (mm) 

dlt  half  of  gear  face  angle  (rad) 

dlta  the  change  in  center  distance  (mm) 

end  the  angle  of  the  meshing  end  (rad) 

f20  the  angle  of  the  first  enveloping  begining  (rad) 

f21  the  angle  of  the  first  enveloping  ended  (rad) 

rl  pitch  radius  of  worm  (mm) 

r2  pitch  radius  of  gear  (mm) 

ral  throat  radius  of  worm  (mm) 

ra2  throat  radius  of  gear  (mm) 

rb  radius  of  base  circle  (mm) 

rfl  root  radius  of  worm  (mm) 

rf2  root  radius  of  gear  (mm) 

reps  the  angle  of  grinding  wheel  turning  to  another  side (rad) 
rgmax  out -radius  of  the  gear  (mm) 

start  the  angle  of  the  meshing  starting  (tad) 

tao  angle  between  two  neighboring  gear  teeth  (tad) 


The  results  which  are  not  called  by  other  subroutines  are  written 
in  the  data  file  under  the  name  of  'dimen.dat'. 


subroutine  dimen ( aO , a2 , ratio , nl , nbao , dl , rfn , bta , 
&  bl,b2,dh, dlta, dlt, f20,f21, 

&  rb , rl , r2 , ral , ra2 , rfl , rf 2 , rgmax , 

&  reps, tao, cmax, cmin, start, end) 

implicit  double  precision(a-h,o-z) 

1  format(2x, 6fl2 .4) 


pitch  diameter  and  radii 

n2=idint(nl*ratio) 
d2=2.0*a2-dl 
rl=dl/2 . 0 
r2=d2/2.0 

circle  pitch,  thickness,  addendum,  dedendum  and  clearance 
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c 


c 


c 

c 

c 


c 


c 


c 

c 

c 


pi— 4 . dO*datan( 1 . dO) 

px=pi*d2/n2 

s2=0.55*px 

sl=px-s2 

xm=px/pi 

gama=datan(px*nl/(pi*dl) ) 
rfx=dtan(rfn)/dcos (gama) 
rfx=datan(rfx) 
drf=px/(2*d2) 

pn=px*dcos (gama) 
ba=.225*pn 
bk=0 . 45*pn 
bt=pn/2 . 0 
bf=bt-ba 
dh=bt-bk 


diameters  and  radius 


db=d2*ds in ( r f x+dr f ) 
dal=dl+2 . 0*ba 
da2=d2+2 . 0*ba 
dfl=dl+2*ba-2*bt 
df2=d2+2*ba-2*bt 
rb=db/2 . 0 
ral=dal/2 . 0 
ra2=da2/2 . 0 
rfl=df 1/2.0 
rf2=df 2/2.0 


rfi=rfx+drf 

tao“2*pi/n2 

psia=. 5*tao*(nbao- .45) 
start=rfi-psia 
end=start+2*psia 
bl=idint(d2*dsin(psia) ) 
b2=idint(.9*dfl) 


rgmax=a2 - dsqrt ( rl*rl -b2*b2/4 .0) 

cmax=a2 - dsqrt ( r f 2*rf 2 -bl*bl/4 . 0) 

cmin=a2 - dsqrt ( ra2*ra2 -bl*bl/4 . 0 ) 

dl ta=cmax+r f 2 - a2 

reps=(2 . 0*rf i+2*sl/d2) 

f20=start-tao 

f21=end+tao 

dl t=das in ( df 1/dl ) 

table  of  the  wormgear  dimensions 


write (1,*) 
write(l,*) 
write (1,1) 
write(l,*) 
write (1,1) 
write (1,*) 
write(l,l) 
write (1,*) 
write (1,1) 
write(l,*) 
write (1,1) 
write(l,*) 


'Dimensions  of  selected  wormgear' 

' aO , a2 , rfn, ratio , db ' 

aO , a2 , rfn, ratio ,  db 

'p2,s2,sl,bk,bt,cc' 

px,s2,sl,bk,bt,dh 

'  dal ,  df  1 ,  da2 ,  df  2 ,  dl ,  d2 ' 

dal,dfl,da2,df2,dl,d2 

' ral , rf 1 , ra2 , rf 2 , rl , r2 ' 

ral , rf 1 , ra2 , rf 2 , rl , r2 

' rf i , rfx , drf , reps , lamda , psia ' 

rf i , rfx , drf  j  reps , gama , psia 

'bta , f20 , f 21 , start , end' 
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write (1,1)  bta,f20,f21,start,end 

write(l,*)  'dlt,bl,b2' 

write(l,l)  dlt,bl,b2 

return 

end 

C 

c  «===============  End  of  Subroutine  DIMEN 

t 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


Subroutine  WORM  called  by  MAIN 

This  subroutine  is  used  to  calculate  the  worm  surface  in  SI,  where 
worm  is  ground  by  a  disc .  To  drive  this  subroutine ,  the  subroutine 
DIMEN  must  be  called  first. 


To  run  this  subroutine  the  input  parameters  are : 

aO  center  distance  of  the  first  enveloping  (mm) 

a2  center  distance  of  wormgear  (mm) 

bl  face  width  of  worm  (mm) 

bt  apex  angle  of  the  generating  plane  (tad) 

cc  clearance  (mm) 

df  step  size  of  fl 

dz  step  size  of  zl 

f20  the  angle  of  the  first  enveloping  begining  (rad) 

f21  the  angle  of  the  second  enveloping  ended  (rad) 

ra2  throat  radius  of  gear  (mm) 

ratio  speed  ratio  of  wormgear 

rb  radius  of  base  circle  (mm) 

rf2  root  radius  of  gear  (mm) 

tb  tan(bt) 


The  variables  in  WORM  are: 

fl  the  running  parameter  in  the  first  enveloping  (rad) 
f2  fl/ratio 

zl  the  coordinate  in  worm  axis  direction 

The  results  which  are  not  called  by  other  subroutines  are  written 
in  the  data  file  under  the  name  of  'worm.dat'. 


subroutine  worm(aO , a2 , ratio , rb ,bt , tb , 

&  f20,f21,bl,ra2,rf2,cc,df ,dz) 

implicit  double  precision(a-h, o-z) 

2  format(2x,6fl0.4) 

Title  of  WORM.dat  file 

write(2,*)  '  The  contact  lines  on  the  worm  surfaces  in  SI' 
write(2,*)  '  . 

write (2,*)  '  fl  |  f2  |  zl  I  I  I 

&  curve 1' 

write(2,*)  '  . - . 

& . ' 

precalculation 


fi=f20*ratio 
fe=f21*ratio 
m=idnint( (fe-f i)/df )+l 
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zi— bl/2.0 
l=idint(bl/dz) 

calculation  of  the  worm  surface  coordinate  in  SI 

do  20  k=l,l 
Z“Z i+(k-l)*dz 
do  10  j“l,m 

fl-fi+(j-l)*df 

f2=fl/ratio 

call  eqw(fl,f2,a0,ratio,rb,tb, 

&  p,q,r,s,pl,ql,rl,sl,p2,q2,r2,s2) 

d,w=p*  (q2*rl  -  ql*r2  )+q*  (pl*r2  -p2*rl)+r*  (p2*ql  -pl*q2  ) 
if(dw.eq.O)  goto  10 

fz=.(p*(q2*sl-ql*s2)+q*(pl*s2-p2*sl)+s*(p2*ql-pl*q2))/dw 
fx=- (q*(r2*sl-s2*rl)+r*(ql*s2-q2*sl)+s*(q2*rl-r2*ql) )/dw 
fy=(p*(r2*sl-s2*rl)+r*(pl*s2-p2*sl)+s*(p2*rl-r2*pl))/dw 

redge=ds  qr  t  (  f x*f x+f3r*f y  ) 
dw=pl*q-p*ql 

x= ( ( ql*r - q*r 1 ) *z - q*sl+ql*s ) /dw 
y=_ ( (pl*r-rl*p)*z-p*sl+pl*s)/dw 

judge  if  the  point  is  in  the  tooth  domain 

if(z.lt.fz)  goto  10 

rw“dsqrt (x*x+y*y ) 

rmax=a2 - dsqr t ( ( r f 2+cc ) **2 - z*z ) 

rmin=a2-dsqrt( (ra2+cc)**2-z*z) 

if(rw.gt.rmax.or.irw.lt.rmin)  goto  10, 

justify  the  tooth  undercut 

if ( redge . gt . rw)  then 

print* , ' the  worm  surface  intersects  with  the  edge  curve ' 
print*,x,y ,z 
endif 

calculation  of  the  worm  surface  curvature 

ff=p2*x+q2*y+r2*z+s2 

rr=r*r 

pel=l . 0+p*p/rr 

pfl=p*q/rr 

pgl=l . 0+q*q/rr 

dw=f f *rr*dsqrt (p*p+q*q+rr ) 

aa==p*rl-r*pl 

bb=q*rl-r*ql 

pll=aa*aa/dw 

pml=aa*bb/dw 

pnl=bb*bb /dw 

curvel=(pel*pnl-2 . 0*pfl*pml+pgl*pll)/(pel*pgl-pfl*pfl) 

write(2,2)  fl,f2,z,x,y,curvel 
10  continue 
20  continue 
return 
end 


End  of  Subroutine  WORM 


c 


C  - ^ - - - - - - - - r-  ■,  r  -  - — . . 

c  Subroutine  EQW  called  by  WORM 
c 

c  This  subroutine  is  used  to  calculate  the  coefficients  of  the  worm 

c  surface  equations . 

C 


c  To  run  this  subroutine  the  input  paraieters  are; 


c  bt  apex  angle  of  the  generating  plane  (rad) 

c  fl  running  parameter  in  the  first  enveloping  (rad) 

c  f2  fl/ratio 

c  rb  radius  of  base  circle  (mm) 

c  tb  tan(bt) 

c 

c  The  output  in  EQW  are : 
c 

c  p,q,r,s  coefficients  of  the  first  equation  in  Eqs(19) 

c  pi,ql,rl,sl  coefficients  of  the  second  equation  in  Eqs(19) 

c  p2,q2,r2,s2  coefficients  of  the  third  equation  in  Eqs(i9) 

c 

c  - 

c 

subroutine  eqw(fl , f2 , aO .ratio , rb , tb, 

&  p,q,r,s,pl,ql,rl,sl,p2,q2,r2,s2) 

implicit  double  precision(a-h,o-z) 
c 

cf2“=dcos(f2) 

sf2=dsin(f2) 

cfl=dcos(fl) 

sfl=dsin(fl) 

c 

c  Here  is  the  equation  (18) 
c 

p=cfl*sf2+tb*sfl 

q=tb*cfl-sfl*sf2 

r=-cf2 

s=rb-a0*sf2 

pl"=tb*cf  l*ratio  -  sf  l*sf  2*ratio+cf  l*cf  2 

ql=-tb*sfl*ratio-sfl*cf2-cfl*sf2*ratio 

rl=sf2 

sl=-a0*cf2 

p2=-tb*sfl*ratio'*ratio-cfl*sf2*ratio*ratio 
&  -cfl*sf2-2*sfl*cf2*ratio 

q2=-tb*cfl*ratio*ratio+sfl*sf2*ratio*ratio 
&  +sfl*sf2-2*cfl*cf2*ratio 

r2-cf2 
s2=a0*sf2 
return 
end 
c 

c . -  -  -END  of  Subroutine  EQW  -  ■  ■  - - - - - - - 

c 

c 

C  - .  - .  -  -  - - 1 - ■  ■■  I  - r.-r'. . . 

c  Subroutine  WORML  called  by  MAIN 
c 

c  This  subroutine  is  used  to  calculate  the  worm  surface  in  SI ,  where 
c  worm  is  ground  by  a  disc.  To  drive  this  subroutine,  the  subroutine 
c  DIMEN  must  be  called  first, 
c 

c  To  run  this  subroutine  the  input  parameters  are: 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


aO  center  distance  of  the  first  enveloping  (mm) 

a2  center  distance  of  wormgear  (mm) 

bl  face  width  of  worm  (mm) 

bt  apex  angle  of  the  generating  plane  (rad) 

cc  clearance  (™o) 

df  step  size  of  fl 

dz  step  size  of  zl 

f20  the  angle  of  the  first  enveloping  begining  (rad) 

f21  the  angle  of  the  second  enveloping  ended  (rad) 

ra2  throat  radius  of  gear  (mm) 

ratio  speed  ratio  of  wormgear 

rb  radius  of  base  circle  (mm) 

reps  angle  of  the  other  side  of  the  worm  surface  (rad) 

rf2  root  radius  of  gear  (mm) 

tb  tan(bt) 


The  variables  in  WORML  are: 

fl  the  running  parameter  in  the  first  enveloping  (rad) 
f2  fl/ratio 

zl  the  coordinate  in  worm  axis  direction 

The  other  results  which  are  not  called  by  other  subroutines  are 
written  in  the  data  file  under  the  name  of  'worml.dat' 


subroutine  worml ( aO , a2 , ratio , rb , b t , tb , 

&  f 20 , f 21 , bl , ra2 , rf 2 , cc , df , dz , reps ) 

implicit  double  precision(a-h,o-z) 

3  format (2x, 6 f 10 .4) 

Title  of  worml.dat  file 

write (3,*)  '  The  contact  lines  on  the  other  side  worm  surfaces' 

write(3,*)  '  - - - - - 7 - 7  7  ' 

write(3,*)  '  fl  I  f2  |  zl  |  xl  ]  ylj 

&  curvel' 

write (3,*)  '  . 

_ ' 


precalculation 

fi=f20*ratio 
fe=f21*ratio 
m=idint( (fe-fi)/df)+l 
zi=-bl/2.0 
l=idint (bl/dz ) 

calculation  of  the  worm  surface  coordinate  in  SI 

do  20  i=l,l 
z=zi+(i-l)*dz 
do  10  j=l,m 
fl=fi+(j-l)*df 
f2=fl/ratio 
f2=f2-reps 

call  eqwl(fl,f2,a0,ratio,rb,tb, 

&  p , q , r , s ,pl , ql , rl , si ,p2 , q2 , r2 , s2) 

dw=p* ( q2*rl - ql*r2 ) +q* (pl*r2 - p2*rl ) +r* (p2*ql - pl*q2 ) 
if  (dw.eq.O)  goto  10 
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fz-- (p*(q2*sl-ql*s2)+q*(pl*s2-p2*sl)+s*(p2*ql-pl*q2) )/dw 
fx»-(q*(r2*sl-s2*rl)+r*(ql*s2-q2*sl)+s*(q2*rl-r2*ql))/dw 
fy“(p*(r2*sl-s2*rl)+r*(pl*s2-p2*sl)+s*(p2*rl-r2*pl))/dw 
redge=dsqr t ( fx*fx+f y*f y ) 


C 

C 


judge  if  the  point  is  in  the  tooth  doiain 


if(z.gt.fz)  goto  10 
dw=pl*q-p*ql 

x=( (ql*r-q*rl)*z-q*sl+ql*s)/dw 
y=- ( (pl*r-rl*p)*z-p*sl+pl*s)/dw 
rw=dsqrt  (x*x+y^) 
rmax=a2-dsqrt( (rf2+cc)**2-z*z) 
nnin=a2-dsqrt( (ra2+cc)**2-z*z) 
if(rw.gt.nnax.or.rw.lt.nnin)  goto  10 
c 

c  judge  the  tooth  undercut 
c 

if(redge.gt.rw)  then 

print*, 'this  point  of  worm  surface  is  on  the  edge  curve' 
print*, x,y,z 
endif 
c 

c  calculation  of  the  worm  surface  curverature 
c 

ff=p2*x+q2*y+r2*z+s2 

rr=r*r 

pel=1.0+p*p/rr 

pfl=p*q/rr 

pgl=l . 0+q*q/rr 

dw=f f*rr*dsqrt (p*p+q*q+rr) 

aa=p*rl-r*pl 

bb=*q*rl-r*ql 

pll=aa*aa/dw 

pml-aa*bb /dw 

pnl=bb*bb /dw 

curvel=(pel*pnl-2 .0*pfl*pml+pgl*pll)/(pel*pgl-pfl*pfl) 
c 

write (3 , 3)  fl , f2+reps ,z ,x,y ,curvel 
10  continue 

20  continue 
return 
end 
c 

c  =======  ■  =  End  of  Subroutine  WORML  ======————=—== 

c 

c 

C  . . . . . ^ - -  - — . . . . . == - - 

c  Subroutine  EQWL  called  by  WORML 
c 

c  This  subroutine  is  used  to  calculate  the  coefficients  of  the  worm 
c  surface . 


c 


c 

To  run 

this  subroutine  the  input  parameters  are: 

c 

aO 

center  distance  of  the  first  enveloping 

(mm) 

c 

bt 

apex  angle  of  the  generating  plane 

(rad) 

c 

fl 

running  parameter  of  the  first  enveloping 

(rad) 

c 

f2 

fl/ratio 

c 

rb 

radius  of  base  circle 

(mm) 

c 

tb 

tan(bt) 
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c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


The  output  in  EQW  are: 


P.qir.s 

pl,ql,rl,sl 

p2,q2,r2,s2 


coefficients 

coefficients 

coefficients 


of 

of 

of 


the  worm 
the  worm 
the  worm 


equation 

equation 

equation 


subroutine  eqwl ( f 1 , f 2 , aO , ratio , rb , tb , 

&  p , q , r , s , pi , ql , rl , si , p2 , q2 , r2 , s2 ) 

implicit  double  precision(a-h,o-z) 

cf2=dcos(f2) 

sf2=dsin(f2) 

cfl=dcos(fl) 

sfl=dsin(fl) 

p=cfl*sf2+tb*sfl 

q=tb*cfl-sfl*sf2 

r=-cf2 

s=-rb-a0*sf2 

pl=tb*cfl*ratio-sfl*sf2*ratio+cfl*cf2 
ql=- tb*sf l*ratio - sf l*cf 2 -cf l*sf 2*ratio 
rl=sf2 
sl=-a0*cf2 

p2=-tb*sfl*ratio*ratio-cfl*sf2*ratio*ratio 
&  -cfl*sf2-2*sfi*cf2*ratio 

q2=-tb*cfl*ratio*ratio+sfl*sf2*ratio*ratio 
&  +sfl*sf2-2*cfl*cf2*ratio 

r2=cf2 
s2=a0*sf2 
return 
end 


End  of  Subroutine  EQWL 


Subroutine  GEAR  called  by  MAIN 

This  subroutine  is  used  to  calculate  the  gear  surface  in  S2,  where 
the  gear  hob  is  ground  by  a  disc.  To  drive  this  subroutine,  the 
subroutine  DIMEN  must  be  called  first. 


To  run  this  subroutine  the  input  parameters  are: 

aO  center  distance  of  the  first  enveloping  (mm) 

a2  center  distance  of  wormgear  (mm) 

bl  face  width  of  worm  (mm) 

b2  face  width  of  gear  (nim) 

bt  apex  angle  of  the  generating  plane  (rad) 

cc  clearance  (mm) 

cmax  maximum  thread  radius  of  the  gear  hob  (mm) 

cmin  minimum  root  radius  of  the  gear  hob  (mm) 

ddd  the  change  in  center  distance  (mm) 

dlt  half  of  gear  face  angle  (rad) 

df  step  size  of  fl 

dz  step  size  of  z2 

f20  the  angle  of  the  first  enveloping  begining  (rad) 

f21  the  angle  of  the  second  enveloping  end  (rad) 

rl  pitch  radius  of  worm  (nim) 

r2  pitch  radius  of  gear  (mm) 

ral  throat  radius  of  worm  (mm) 
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c 

ra2 

throat  radius  of  gear 

(mm) 

c 

ratio 

speed  ratio  of  wormgear 

c 

ratio2 

speed  ratio  of  wormgear 

c 

rb 

radius  of  base  circle 

(mm) 

c 

rfl 

root  radius  of  worm 

(mm) 

C 

rf2 

root  radius  of  gear 

(mm) 

C 

niiax2 

out -radius  of  the  gear 

(min) 

running  angle  of  the  second  enveloping  start 
fiiMiiii  Msk  jf  fk  asaak  anvalanliis  auk/ 


(rad) 

ha# 


til  fiiMidi  Mjia  at  tha  aaaana  aavaiaaiaa  aiiflafl  IM) 

tb  tan(bt) 

The  variables  in  GEAR  are : 

tl  the  running  parameter  in  the  second  enveloping  (rad) 
t2  tl/ratio 

z2  the  coordinate  in  gear  axis  direction 

The  results  which  are  not  called  by  other  subroutines  are  written 
in  the  data  file  under  the  name  of  'gear.dat'. 


subroutine  gear(aO ,a2 , ratio ,ratio2 , rb , tb ,bt , rmax2 , cmax, cmin, 

&  ra2 , rf 2 , r2 , ral , rf 1 , rl , cc , ddd ,bl ,b2 , 

&  f20,f21,df,t20,t21,dt,dlt,dz) 

implicit  double  precision(a-h,o-z) 

4  format(2x, 6fl2 .4) 

Title  of  Gear.dat  file 

write(4,*)'  This  is  the  point  coordinates  on  gear  surface; (  when 
&  tl=fl,  the  point  is  on  the  first  contact  line)' 

write (4,*)  ' . 

& . ' 

write(4.*)  '  tl  [  fl  |  x2  |  y2  | 

&  z2  I  r' 

write (4,*)  '  . - . 

& . ' 

l=idint(b2/dz)+l 
m=idint((t21-t20)/dt)+l 
df 1-df /lO . 0 

n=idint(4. dO*dlt/dfl)+l 

precalculation 

do  10  i=l,m 

t2=t20+(i-l)*dt 

tl=t2*ratio2 

calculation  of  the  curved  part  of  the  gear  surface  in  S2 
do  20  j=l,n 

fl=(tl-2 . 0*dlt)+(i -l)*dfl 
f2=fl/ratio 

if(f2.1t.f20.or.f2.gt.f21)  goto  20 

call  second( t2 , f 2 , aO , a2 , ratio , ratio2 , rb , tb ,bt , 

&  xl,yl,zl,x2,y2,z2,xn,yn,zn, 

&  curve 1 , curve 2 , curve , dw ) 


judge  if  the  point  is  in  the  tooth  domain 
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if(dw.eq.O)  goto  20 

if (curve. It .0)  goto  20 

if(z2.gt.b2/2.or.z2.1t.-b2/2)  goto  20 

r g=ds  qr t ( x2*x2+y 2*y 2 ) 

rmin= ( a2+ddd) - dsqrt ( cmax**2 - z2*z2 ) 

nnax=(a2+ddd)  -  dsqrt  ( cinin**2  -  z2*z2  ) 

if  ( rmax . gt . rmax2 )  then 

rmax=rmax2 

endif 

if(rg.gt.rmax.or.rg.lt.rmin)  goto  20 

write (4 , 4) tl , f 1 , x2 , y2 , z2 , rg 
20  continue 


c 

c 

c 

c 


c 

c 

c 


c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


if (a0.ne.a2.or.ratio.ne.ratio2)  goto  10 

calculation  of  the  planar  part  of  the  gear  surface  in  S2 
(  for  the  case  of  fl=tl  in  the  conventional  form  ) 

do  30  k=l,l 

z2=-b2/2+(k-l)*dz 

fl=tl 

f2=t2 

if(f2.1t.f20.or.f2.gt.f21)  goto  30 

call  first(t2,f2,a0,a2, ratio , rat io2 , rb , tb , b t , z 2 , 

&  xl,yl,zl,y2,x2,xn,yn,zn, 

&  curve 1 , curve 2 , curve , dw) 

judge  if  the  point  is  in  the  tooth  domain 

if(dw.eq.O)  goto  30 

if (curve, It. 0)  goto  30 

rg=ds  qr t (x2*x2+y 2*y 2 ) 

rmin=  ( a2+ddd) - ds qr t ( cmax**2 - z2*z 2 ) 

rmax= ( a2+ddd) - dsqrt ( cmin**2 - z2*z2 ) 

if  ( rmax . gt . rmax2 )  then 

rmax==rmax2 

endif 

if(rg,gt.rmax.or.rg.lt,rmin)  goto  30 

write (4 , 4) tl , f 1 , x2 , y2 , z2 , rg 
30  continue 

10  continue 
return 
end 


End  of  Subroutine  GEAR 


Subroutine  SECOND  called  by  Subroutine  GEAR 

This  subroutine  is  used  to  calculate  the  coefficiences  of  the  gear 
surface  equations  Eq(28)  in  techanical  report. 


To  run  this  subroutine  the  input  parameters  are: 

aO  center  distance  of  the  first  enveloping  (mm) 
a2  center  distance  of  wormgear  (mm) 
btr  apex  angle  of  the  generating  plane  (rad) 
f2  running  parameter  in  the  first  enveloping  (rad) 
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c 

c 

c 

c 

C 

C 

c 

{ 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 


c 

c 

c 


c 

c 

c 


ratio  speed  ratio  of  wormgear 
ratio2  speed  ratio  of  wormgear 

rb  radius  of  base  circle  (mm) 
t2  running  parameter  in  the  second  enveloping  (rad) 
tb  tan(bt) 


curve 1 
curve 2 
curve 
dw 

xl ,yl , zl 
x2,y2,z2 
xn , yn , zn 


worm  surface  curvature  at 
gear  surface  curvature  at 
relative  curvature  at  the 
index  (  dw=0  the  point  is 
contact  point  coordinates 
contact  point  coordinates 
contact  point  coordinates 


the  contact  point 
the  contact  point 
contact  point 
a  singular  point  ) 
in  SI 
in  S2 
in  Sn 


subroutine  second(t2 , f2 ,a0 ,a2 .ratio ,ratio2 , rb , tb ,btr , 

&  xl,yl,zl,x2,y2,z2,xn,3nri,zn, 

&  curve 1 , curve  2 , curve , dw ) 

implicit  double  precision(a-h,o-z) 

calculate  the  coordinates  of  the  point  for  given  fl,  tl 

tl=t2*ratio2 

fl=f2*ratio 

call  eqg( tl , t2 , f 1 , f 2 , aO , a2 , ratio , ratio2 , tb , rb , btr , 

&  p , q , r , s , pi , ql , rl , si , p2 , q2 , r2 , s2 ) 

dw=p*(q2*rl-qi*r2)+q*(pl*r2-p2*rl)+r*(p2*ql-pl*q2) 
if  (dw.eq.O)  goto  10 

x2 — (q*(r2*sl-s2*rl)+r*(ql*s2-q2*sl)+s*(q2*rl-r2*ql))/dw 
y2=(p*(r2*sl-s2*rl)+r*(pl*s2-p2*sl)+s*(p2*rl-r2*pl) )/dw 
z2=- (p*(q2*sl-ql*s2)+q*(pl*s2-p2*sl)+s*(p2*ql-pl*q2) )/dw 

ctl=dcos (tl) 

stl=dsin(tl) 

ct2=dcos (t2) 

st2=dsin(t2) 

cf2=dcos(f2) 

sf2=ds in(f2) 

cfl=dcos(fl) 

sfl=dsin(fl) 

sft=dsin(fl-tl) 

cft“dcos(fl-tl) 

Coordinate  transformation 

xl=-x2*ctl*ct2+y2*ctl*st2-z2*stl+a2*ctl 

yl=x2*stl*ct2-y2*stl*st2-z2*ctl-a2*stl 

zl=-x2*st2-y2*ct2 

xn=x2*ct2-y2*st2 

yn=x2*st2+y2*ct2 

zn=z2 

curvature  of  the  gear  surface 


rr=r*r 

pe2=l+p*p/rr 

pf2=p*q/rr 

pg2=l+q*q/rr 
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d=q*rl-r*ql 
c=p*r2-r*p2 
e=q*r2-r*q2 


rto=ratio*ratio 

rto2=ratio2*ratio2 

f  fx>=sf  2*ct2*cf  t+tb*ct2*sf  t+2 . 0*cf  2*ct2*sf  t/ratio 
&  +(sf2*ct2*cft-cf2*st2)/rto 

f fy-sf 2*st2*cf t+tb*st2*sf t+2 . 0*cf2*st2*sf t/ratio 
&  +(sf2*st2*cft+cf2*ct2)/rto 

f f z=s  f 2*s  ft - tb*cf t - 2 . 0*cf 2*cf t/rat io+sf 2*s  f t/rto 
ffd=a2*(sf2*cft+tb*sft+2 . 0*cf2*sft/ratio+sf2*cft/rto) 

&  -a0*sf2/rto 

f f=f fx*x2 - f fy*y2 - f f z*z2 - f f d 

ftx=sf2*st2*sft-tb*st2*cft+ratio2*(sf2*ct2*cft+tb*ct2*sft) 

&  +(sf2*ct2-cf2*st2*cft+ratio2*cf2*ct2*sft)/ratio 

fty=-sf2*ct2*sft+tb*ct2*cft+ratio2*(sf2*st2*cft+tb*st2*sft) 
&  +(sf2*st2+cf2*ct2*cft+ratio2*cf2*st2*sft)/ratio 
f tz=ratio2* ( s  f 2*sf t - tb*cf t - cf2*cf t/ratio) 
ftd=ratio2*a2*(sf2*cft+tb*sft+cf2*sft/ratio) 
ft=-ftx*x2+fty^2+ftz*z2+ftd 

ttx=-p+2 . 0*ratio2*(sf2*st2*sf t-tb*st2*cft) 

&  +rto2*(sf2*ct2*cft+tb*ct2*sft) 

tty=-q+2.0*ratio2*(sf2*ct2*sft-tb*ct2*cft) 

&  -rto2*(sf2*st2*cft+tb*st2*sft) 

ttz=rto2*(tb*cft-sf2*sft) 
ttd=rto2*a2*(sf2*cft+tb*sft) 
tt=ttx*x2+tty*y2+ttz*z2 - ttd 

g=ff*tt-ft*ft 
cm=dcos (btr ) / ( g*rr ) 
pl2=cni*(b*b*tt-2*b*c*ft+c*c*ff) 
pm2=cm* (b*d*tt- (b*e+c*d) *f t+c*e*f f ) 
pn2=cin*  (  d*d*tt  -  2'*d*e*f  t+e*e*f  f  ) 

Curve2=(pe2*pn2-2*pf2*pm2+pg2*pl2)/(pe2*pg2-pf2*pf2) 

curvature  of  the  worm  surface 

pel“l . 0+p*p/rr 

pfl=p*q/rr 

pgl=l . 0+q*q/rr 

dw=f  f *r r*ds  qr t (p*p+q*q+rr ) 

aa=p*rl-r*pl 

bb=q*rl-r*ql 

pll=aa*aa/dw 

pml=aa*bb/dw 

pnl=bb*bb/dw 

curvel=(pel*pnl-2 . 0*pfl*pml+pgl*pll)/(pel*pgl-pfl*pfl) 
relative  curvature 

cuirve=curve2  -  curve  1 


10  return 
end 


End  of  Subroutine  SECOND 


c 


Subroutine  EQG  called  by  Subroutine  FIRST  and  SECOND 

This  subroutine  is  used  to  calculate  the  coefficients  of  the  gear 

surface  equations. 

lo  ri  this  subroutine  the  input  paraieters  are; 


itef  itshnce  of  t:L  flufit 


aO  center  distance  or  tke  riHfifc 

a2  center  distance  of  the  first  enveloping 

bt  apex  angle  of  the  generating  plane 
fl  running  parameter  of  the  first  enveloping 
f2  fl/ratio 

tl  running  parameter  of  the  second  enveloping 
t2  tl/ratio 

ratio  speed  ratio  of  wormgear 
ratio2  speed  ratio  of  wormgear 
rb  radius  of  base  circle 
tb  tan(bt) 

The  output  in  EQW  are : 

p,q,r,s  coefficients  in  Eqs(29) 

pi,ql,rl,sl  coefficients  in  Eqs(30) 
p2,q2,r2,s2  coefficients  in  Eqs(31) 


(mm) 

(rad) 

(rad) 

(rad) 


subroutine  eqg( tl , t2 , f 1 , f 2 , aO , a2 , ratio , ratio2 , tb , rb ,btr , 

&  p , q , r , s ,pl , ql , rl , si , p2 , q2 , r2 , s2 ) 

implicit  double  precision(a-h,o-z) 

ctl— dcos(tl) 

stl=dsin(tl) 

ct2=dcos(t2) 

st2=dsin(t2) 

cf2=dcos(f2) 

sf2=dsin(f2) 

cfl=dcos(fl) 

sfl=dsin(fl) 

sft=dsin(fl-tl) 

cft=dcos(fl-tl) 

Here  is  the  Eqs(29)-(30) 

p=-sf2*ct2*cft-tb*ct2*sft+cf2*st2 
q=sf2*st2*cft+tb*st2*sft+cf2*ct2 
r=sf2*sft-tb*cft 
s=a2*(sf2*cft+tb*sft) -a0*sf2+rb 

pl=sf2*ct2*sft-tb*ct2*cft- (sf2*st2+cf2*ct2*cft)/ratio 
ql=-sf2*st2*sft+tb*st2*cft- (sf2*ct2-cf2*st2*cft)/ratio 
rl=sf2*cft+tb*sft+cf2*sft/ratio 

sl=-a2*(sf2*sft- tb*cft-cf2*cft/ratio) -a0*cf2/ratio 

p2=q-ratio2*(sf2*ct2*sft-tb*ct2*cft) 

q2--p+ratio2*(sf2*st2*sft-tb*st2*cft) 

r2=-ratio2*(sf2*cft+tb*sft) 

s2=a2*ratio2*(sf2*sft-tb*cft) 

return 

end 


End  of  Subroutine  EQG 
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c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 


c 

c 
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Subroutine  FIRST  called  by  GEAR 

This  subroutine  is  used  to  calculate  the  coefficients  of  the  gear 
surface  equations  and  curvatures. 


To  run  this  subroutine  the  input  parameters  are: 

aO  center  distance  of  the  first  enveloping  (mm) 

a2  center  distance  of  wormgear  (mm) 

btr  apex  angle  of  the  generating  plane  (rad) 

f2  running  parameter  in  the  first  enveloping  (rad) 

ratio  speed  ratio  of  wormgear 
ratio2  speed  ratio  of  wormgear 

rb  radius  of  base  circle  _  (mm) 

t2  running  parameter  in  the  second  enveloping  (rad) 

tb  tan(bt) 

z2  coordinate  of  the  point  in  S2 


The  output  in  FIRST  are: 


curve 1 
curve2 
curve 
dw 

xl,yl,zl 
x2,y2 
xn ,  3m , zn 


worm  surface  curvature  at 
gear  surface  curvature  at 
relative  curvature  at  the 
index  (  dw=0  the  point  is 
contact  point  coordinates 
contact  point  coordinates 
contact  point  coordinates 


the  contact  point 
the  contact  point 
contact  point 
a  singularity  point 
in  SI 
in  S2 
in  Sn 


) 


subroutine  f irs t ( t2 , f 2 , aO , a2 , ratio , ratio2 , rb , tb , btr , z2 , 
&  xl,yl,zl,y2,x2,xn,3m,zn, 

&  curve! , curve2 , curve , dw) 

implicit  double  precision(a-h,o-z) 

calculation  of  the  coordinates  of  the  point 
at  the  given  position  fl,  tl 

tl=t2*ratio2 

fl=f2*ratio 

call  eqg ( tl , t2 , f 1 , f 2 , aO , a2 , ratio , ratio2 , tb , rb , b tr , 

&  p , q , r , s , pi , ql , rl , si ,p2 , q2 , r2 , s2) 

dw=pl*q-p*ql 
if  (dw.eq.O)  goto  10 
x2= ( ( ql^r - q*r 1 ) *z2 - q*s l+ql*s ) /dw 
y2=- ( (pl*r-p*rl)*z2-p*sl+pl*s)/dw 

ctl=dcos(tl) 

stl=dsin(tl) 

ct2=dcos(t2) 

st2=dsin(t2) 

cf2=dcos(f2) 

sf2=dsin(f2) 

cfl=dcos(fl) 

sfl=dsin(fl) 

sft=dsin(fl-tl) 

cft=dcos(fl-tl) 

Coordinate  transformation 
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xl“-x2*ctl*ct2+y2*ctl*st2-z2*stl+a2*ctl 

yl=x2*stl*ct2-y2*stl*st2-z2*ctl-a2*stl 

zl'=-x2*st2-y2*ct2 

xn=x2*ct2 -y2*st2 

yn=x2*st2+y2*ct2 

zn=z2 


of  ttie  Eoar  surface 


rr=r*r 

pe2=l+p*p/rr 

pf2=p*q/rr 

pg2=l+q*q/rr 

b-p*rl-r*pl 

d-=q*rl-r*ql 

c-b 

e=d 


rto=ratio*ratio 

rto2=ratio2*ratio2 

f fx=sf 2*ct2*cf t+tb*ct2*sf t+2 . 0*cf 2*ct2*sf t/ratio 
&  +(sf2*ct2*cft-cf2*st2)/rto 

f fy=sf 2*st2*cf t+tb*st2*sf t+2 . 0*cf 2*st2*sf t/ratio 
&  +(sf2*st2*cft+cf2*ct2)/rto 

ffz=sf2*sft-tb*cft-2.0*cf2*cft/ratio+sf2*sft/rto 

ffd=a2*(sf2*cft+tb*sft+2.0*cf2*sft/ratio+sf2*cft/rto) 

&  -aO*sf2/rto 

f f=f fx*x2 - f fy*y2 - f f z*z2 - f f d 

ftx=sf2*st2*sft- tb*st2*cft+ratio2*(sf2*ct2*cft+tb*ct2*sft) 

&  +(sf2*ct2-cf2*st2*cft+ratio2*cf2*ct2*sft) /ratio 

fty=-sf2*ct2*sft+tb*ct2*cft+ratio2*(sf2*st2*cft+tb*st2*sft) 
&  +(sf2*st2+cf2*ct2*cft+ratio2*cf2*st2*sft) /ratio 
f tz=ratio2* ( sf 2*sf t - tb*cf t- cf2*cf t/ratio) 
ftd=ratio2*a2*(sf2*cft+tb*sft+cf2*sf t/ratio) 
ft=-ftx*x2+fty*y2+ftz*z2+ftd 

ttx=-p+2 . O*ratio2*(sf2*st2*sft- tb*st2*cft) 

&  +rto2*(sf2*ct2*cft+tb*ct2*sft) 

tty=-q+2 .O*ratio2*(sf2*ct2*sft-tb*ct2*cft) 

&  -rto2*(sf2*st2*cft+tb*st2*sft) 

ttz=rto2*(tb*cft-sf2*sft) 
ttd=rto2*a2*(sf2*cft+tb*sft) 
tt=ttx*x2+tty*y2+ttz*z2 - ttd 

g=ff*tt-ft*ft 

cm=dcos (btr)/(g*rr) 

pl2=0 

pm2=0 

pn2=0 

Curve2=(pe2*pn2-2*pf2*pm2+pg2*pl2)/(pe2*pg2-pf2*pf2) 

curvature  of  the  worm  surface 

pel=l . 0+p*p/rr 

pfl=p*q/rr 

pgl=l . 0+q*q/rr 

dw=f f*rr*dsqrt (p*p+q*q+rr) 

aa=p*rl-r*pl 

bb=q*rl-r*ql 
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pll=aa*aa/dw 

pinl=aa*bb/dw 

pnl=bb*bb/dw  ,  . 

curvel=(pel*pnl-2 . 0*pfl*pml+pgl*pll)/(pel*pgl-pfl*pfl) 

relative  curvature 

curve=curve2 - curvel 

D  return 
end 

■END  of  Subroutine  FIRST 


Subroutine  CONTACT  called  by  MAIN 

This  subroutine  is  used  to  calculate  the  real  contact  lines  in  Sn 
for  the  given  position  tn2.  To  drive  this  subroutine,  the 
subroutine  DIMEN  must  be  called  first  and  the  given  position  must 
be  given  previously. 


To  run  this  subroutine  the  input  parameters  are: 

aO  center  distance  of  the  first  enveloping 

a2  center  distance  of  wormgear 

bl  face  width  of  worm 

b2  face  width  of  gear 

btr  apex  angle  of  the  generating  plane 

dh  clearance 

cmax  maximum  thread  radius  of  the  gear  hob 

cmin  minimum  root  radius  of  the  gear  hob 

ddd  the  change  in  center  distance 

dlt  half  of  gear  face  angle 

df  step  size  of  fl 

dz  step  size  of  z2 

end  the  angle  of  the  meshing  ended 

f20  the  angle  of  the  first  enveloping  begining 

f21  the  angle  of  the  second  enveloping  ended 

np  number  of  meshing  teeth 

rl  pitch  radius  of  worm 

r2  pitch  radius  of  gear 

ral  throat  radius  of  worm 

ra2  throat  radius  of  gear 

ratio  speed  ratio  of  wormgear 

ratio2  speed  ratio  of  wormgear 

rb  radius  of  base  circle 

rfl  root  radius  of  worm 

rf2  root  radius  of  gear 

rmax2  outradius  of  the  gear 

start  the  angle  of  the  meshing  starting 

t20  running  angle  of  the  second  enveloping  start 

t21  running  angle  of  the  second  enveloping  ended 

tb  tan(bt) 

tn2  the  given  position 

tt2  angle  between  two  neighboring  gear  teeth 


(mm) 

(mm) 

(mm) 

(mm) 

(rad) 

(mm) 

(mm) 

(mm) 

(mm) 

(rad) 


(rad) 

(rad) 

(rad) 


(mm) 

(mm) 

(mm) 

(mm) 

(rad) 

(rad) 

(rad) 


(rad) 


The  output  of  the  subroutine  are: 

contc  the  length  of  contact  lines  of  each  meshing  tooth 
nc  the  niimber  of  real  meshing  pair 
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The  variables  in  CONTACT  are: 

tl  the  running  parameter  in  the  second  enveloping 

t2  tl/ratio 

z2  the  coordinate  in  gear  axis  direction 


(rad) 


tke  common  excliange  iimensions  are! 

x,y,z  save  the  coordinates  of  contact  point  on  meshing  pair 
nn  maximum  dimension  number  of  the  exchange  dimension 

rg  save  the  radius  of  contact  point  to  gear  axis 

The  results  which  are  not  called  by  other  subroutines  are  written 
in  the  data  file  under  the  name  of  'contc.dat'. 

subroutine  contact (aO , a2 , ratio , ratio2 , np , rb , tb , btr , dh ,  ddd , 

&  rl , ral , rf 1 , r2 , ra2 , rf 2 , rmax2 , cmax , cmin , bl , b2 , 

&  f 20 , f 2i , df , t20 , t2i , tt2 , tn2 , start , end , dlt , dz , 

&  contc,nc,nn,x,y,z,rg) 

implicit  double  precision(a-h,o-z) 
dimension  x(nn) ,y(nn) ,z(nn) ,rg(nn) ,contc(np) 

5  format(2x,6fl0.4) 

Title  of  contc.dat  file 

write(5,*)'  The  length  of  contact  lines  on  each  meshing  teeth' 

write(5,*)  '  . - . - . ' 

write (5,*)  '  tnl  |  first  ]  second  |  total  ' 

write(5,*)  '  . ' 

precalculation 

pi=4 . d0*datan(l . dO) 
df l=df/100 . 0 
dz2=dz/10 . 0 
nc=0 

l=idint (b2/dz2 )+l 
m=idint(f21/tt2)+l 
n=idint(4 . 0*dlt/dfl)+l 


calculation  of  the  coordinates  of  the  contact  point 

do  10  i=l,m 

t2=tn2+(i-l)*tt2 

if (t2.lt. start. or. t2.gt. end)  goto  10 

tl'=ratio2*t2 

ctl-=dcos  (tl) 

stl=dsin(tl) 

ct2=dcos(t2) 

st2=dsin(t2) 

ind=0 

do  20  j=l,n 

f 1= ( tl -  2 . 0*dl t ) + ( j - 1 ) *df 1 
f2=fl/ratio 

if (f2 . It . start . or . f2 . gt . end)  goto  20 
call  second( t2 , f 2 , aO , a2 , ratio , ratio2 , rb , tb ,bt , 
&  xl,yl,zl,x2,y2,z2,xn,yn,zn, 
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c 
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&  curvel,curve2, curve, dw) 

judge  if  the  point  is  in  the  tooth  domain 

if(dw.eq.O)  goto  20 
if (curve. It. 0)  goto  20 

if(z2.1t.-b2/2.0.or.z2.gt.b2/2.0)  goto  20 
rg2=dsqrt (x2*x2+y2*y2 ) 

rn2=(a2+ddd)  -dsqrt(  (cmax-dh)'**2-z2*z2) 
rm2“ ( a2+ddd) - dsqrt ( cmin**2 - z  2*z2 ) 
if  (rm2.gt.rmax2)  then 
rm2=rmax2 
endif 

if (rg2. gt.rm2.or.rg2.lt. m2)  goto  20 

if(zl.gt.bl/2.0.or.zl.lt.-bl/2.0)  goto  20 
rwl=dsqrt (xl*xl+yl*yl ) 
rnl=a2 - dsqrt ( ( ra2+dh) **2 - zl*zl) 
rml=a2- dsqrt ( (rf2+dh)**2-zl*zl) 
if(rwl.gt.rml.or.rwl.lt.rnl)  goto  20 
c 

ind=ind+l 

if (ind. gt .nn)  then 

print*,'  demensions  nn  is  not  big  enough  for  contact' 
print*,'  increase  nn  or  df  in  main  program.  ' 
print* , ind 
pause 
endif 
x(ind)=xn 
y(ind)’=yn 
z(ind)=z2 
rg(ind)=rg2 
20  continue 
c 

c  calculation  of  the  length  of  the  second  contact  line 
c 

if  (ind.ne.O)  then 

rcl2=dsqrt( (x(ind) -x(l) )**2+(y(ind) -y(l) )**2 
&  +(z(ind)-z(l))**2) 

endif 

do  30  j=l,ind 
x(j)=0 

y(j)=o 

z(j)=0 

30  continue 
c 

if(a0.ne.a2. or. ratio. ne.ratio2)  goto  10 
c 

c  calculation  of  the  coordinates  of  the  contact  point 
c  (this  is  for  the  situation  where  fl=tl) 
c 

ind=0 

do  60  k=l,l 

z2=-b2/2.0+(k-l)*dz2 

fl=tl 

f2=t2 

if(f2.1t.f20.or.f2.gt.f21)  goto  60 
c 

call  f irst ( t2 , f 2 , aO , a2 , ratio , ratio2 , rb , tb ,bt , z2 , 
xl,yl,zl,y2,x2,xn,yn,zn, 
curvel , curve 2 , curve , dw) 


& 

6c 
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c 

c 

c 


C 


c 


c 


c 

c 

c 


c 

c 

c 


c 


c 

c 


c 

c 


judge  if  the  point  is  in  the  tooth  domain 

if(dw.eq.O)  goto  60 

if (curve. It. 0)  goto  60 

if(zl.gt.bl/2.0.or.zl.lt.-bl/2.0)  goto  60 

rwl-dsqrt (xl*xl+yl*yl ) 
rml-a2-dsqrt( (rf2+dh)*(rf2+dh) -zl*zl) 
ml-a2-dsqrt( (ra2+dh)*(ra2+dh) -zl*zl) 
if(rwl.gt.rml.or.rwl.lt.rnl)  goto  60 

if(z2.1t. -b2/2.0.or.z2.gt.b2/2.0)  goto  60 

rg2»=dsqrt(x2*x2+y2*y2) 

rn2=(a0+ddd) -dsqrt( (cmax-dh)**2-z2*z2) 

rm2= ( aO+ddd ) - ds  qr t ( cmin**2 - z  2*z  2 ) 

if  (rm2.gt.rmax2)  then 

rm2=rmax2 

endif 

if(rg2.gt.rm2.or.rg2.1t.rn2)  goto  60 
ind=ind+l 

if(ind.gt.nn)  then 

print*,'  demensions  nn  is  not  big  enough  for  contact' 
print*,'  increase  nn  or  dz  in  main  program.  ' 
print* , ind 
pause 
endif 
x(ind)=xn 
y(ind)=yn 
z(ind)=z2 
rg(ind)=rg2 
60  continue 

calculation  of  the  length  of  the  second  contact  line 
if  (ind.ne.O)  then 

rcll=dsqrt( (x(ind) -x(l) )**2+(y(ind) -y(l) )**2 
&  +(z(ind)-z(l))**2) 

endif 


calculation  of  the  sum  of  two  contact  lines  on  the  meshing  tooth 

if(rcll.eq.0.and.rcl2.eq.0)  goto  10 
nc=nc+l 

contc(nc)=rcll+rcl2 

write(5,5)  tl , rcll , rcl2 ,contc(nc) 

rcll=0 

rcl2=0 


70 


do  70  j=l,ind 
x(j)=0 

y(j)=o 

z(j)=0 

continue 


10  continue 


return 

end 


END  of  Subroutine  CONTACT 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Subroutine  PITCH  called  by  MAIN 

This  subroutine  is  used  to  calculate  the  coordinates  of  the  worm 
pitch  points,  its  normal  direction  and  relative  velocity  at  this 
point  in  Sj .  All  the  calculations  are  for  the  given  position  tn2 
To  drive  this  subroutine  DIMEN  must  be  called  first  and  the  given 
position  must  be  input  previously. 


To  run  this  subroutine  the  input  parameters  are; 

aO  center  distance  of  the  first  enveloping  (mm) 

a2  center  distance  of  wormgear  (mm) 

bl  face  width  of  worm  (mm) 

b2  face  width  of  gear  (mm) 

btr  apex  angle  of  the  generating  plane  (rad) 

dh  clearance  (JJ™) 

cmax  maximum  throat  radius  of  the  gear  hob  (mm) 

cmin  minimum  root  radius  of  the  gear  hob  (mm) 

ddd  the  change  in  center  distance  (mm) 

dlt  half  of  gear  face  angle  (rad) 

df  step  size  of  fl 

dz  step  size  of  z2 

end  the  angle  of  the  meshing  ended  (rad) 

ep.epp  convergence  criterion  for  root  searching 
f20  the  angle  of  the  first  enveloping  begining  (rad) 

f21  the  angle  of  the  second  enveloping  ended  (rad) 

np  nvimber  of  meshing  teeth 

rl  pitch  radius  of  worm  (mm) 

r2  pitch  radius  of  gear  (mm) 

ral  throat  radius  of  worm  (mm) 

ra2  throat  radius  of  gear  (mm) 

ratio  speed  ratio  of  wormgear 
ratio2  speed  ratio  of  wormgear 

rb  radius  of  base  circle  (mm) 

rfl  root  radius  of  worm  (mm) 

rf2  root  radius  of  gear  (mm) 

rmax2  outradius  of  the  gear  (mm) 

start  the  angle  of  the  meshing  starting  (rad) 


t20  running  angle  of  the  second  enveloping  start  (rad) 

t21  running  angle  of  the  second  enveloping  ended  (rad) 

tb  tan(bt) 

tn2  the  given  position 

tt2  angle  between  two  neighboring  gear  teeth  (rad) 


The  output  of  the  subroutine  are: 


nc 

tp2 

fp2 

pxj ,pzj 
xpn , ypn , zpn 

vx,vy ,vz 


the  number  of  real  meshing  teeth 
the  position  of  each  meshing  pair 
the  f2  value  of  pitch  points 
the  coordinates  of  pitch  points  in  Sj 
the  normal  direction  components  of 
pitch  point  in  Sj 

the  relative  velocity  components  of 
pitch  point  in  Sj 


The  variables  in  PITCH  are: 

tl  the  running  parameter  in  the  second  enveloping  (rad) 
t2  tl/ratio 
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c  zl  the  coordinate  in  worm  axis  direction 

c 

c 

c  The  common  exchange  dimensions  are; 

c  ff  save  the  fl  values  of  the  worm  pitch  points 

c  xin,zm  save  the  coordinates  of  the  worm  mid-plane  profile 

c  rr  save  the  radius  of  worm  pitch  points 

c  nn  laxiii  diiension  Biber  of  the  exchange  diiension 

c 

c 

c  The  results  which  are  not  called  by  other  subroutines  are  written 
c  in  the  data  file  under  the  name  of  'pitch.dat'. 

c - - - - - . . . 

c 

subroutine  pitch(aO ,a2 , ratio, ratio2 ,np ,rb , tb ,btr , dh, ddd, 

&  rl , ral , rf 1 , r2 , ra2 , rf 2 , rmax2 , cmax , cmin , bl , b2 , 

&  f 20 , f 21 , df , t20 , t21 , tt2 , tn2 , start , end , dlt , dz , 

&  nn,xm,zm,ff ,rr,ep,epp, 

&  tp2,fp2,pxj ,pzj ,xpn,ypn,zpn,vx,vy,vz,nc) 

implicit  double  precision(a-h, o-z) 

dimension  xm(nn) ,zm(nn) ,ff(nn) ,rr(nn) ,tp2(np) ,fp2(np) , 

&  pxj (np) ,pzj (np) ,xpn(np) ,ypn(np) ,zpn(np) , 

&  vx(np) ,vy(np) ,vz(np) 

6  format(2x, 6fl2 .4) 
c 
c 

c  Title  of  contc.dat  file 
c 

write(7,*)  '  pitch  point  and  normal  direction  ' 

write(7,*)  '  . - . - . .  ' 

write(7,*)  '  tnl  ]  xj  |  zj  |  xpn  | 

&  ypn  zpn' 

write(7,*)  '  . 

& . ' 

write (8,*)  '  pitch  point  and  relative  velocity  direction' 

write (8,*)  '  . - . - . - . -- 

write (8,*)  '  tnl  |  xj  |  zj  |  rvx  | 

&  rvy  1  rvz' 

write  (8,*)  '  - - - . . 

& . ' 

c 

c  Precalculation 
c 

nc=0 
yl=0 . 0 

pi-4 . dO*datan(l . dO) 
m=idint ( f 21/tt2 ) +1 
startl=start*ratio 
endl=end*ratio 
dfl=df/100.0 

ml=idint(pi/(2 . 0*dfl) )+l 
c 

c  Worm  profile  on  the  mid-plane  (yj=0) 
c 

do  50  k=l,m 
ind=0 

f n2=tn2+ (k - 1 ) *tt2 
fnl=fn2*ratio2 

if(fnl.lt.startl.or.fnl.gt.endl)  goto  50 
c 

do  10  i=l,ml 

fl=fnl-pi/4.0+(i-l)*dfl 


124 


call  wormt ( f 1 , f nl , y 1 , aO , ratio , rb , tb , b tr , xl , zl , mw) 
if(niw.eq.l)  goto  10 
if(xl.lt.O)  goto  10 
if (dabs(zl) .gt.rb)  goto  10 
c 

rwl=dsqrt (xl*xl+yl*yl ) 
ml=aO-dsqrt ( (ra2+dh)*(ra2+dh) -zl*zl) 
nnl=a0-dsqrt((rf2+dh)*(rf2+dh)-zl*zl) 
if (rwl. gt.rml.or.rwl.lt. ml)  goto  10 
c 

ind=ind+l 

if ( ind . gt . nn)  then 

print*, 'the  dimension  is  not  big  enough  in  PITCH' 

print*, ' increase  nn  in  main' 

print*, ind 

pause 

endif 

xm(ind)=xl 

zm(ind)=zl 

rr ( ind) =dsqr t ( ( aO -xl) * ( aO -xl)+zl*zl) 
ff(ind)=fl 
10  continue 
c 


c 

c 

c 


c 


Find  worm  pitch  points 
jt=l 

dr=dabs(rr(l) -r2) 
do  20  j=2,ind 

if(dabs(rr(j)-r2) .It.dr)  then 

dr=dabs (rr ( j ) - r2 ) 

jt=i 

endif 

20  continue 

if(ff(jt) .gt.endl.or.ff(jt) .It.startl)  goto  25 


fl=ff(jt-l) 

fr=ff(jt+l)  ,  ^ 

call  roo t ( f 1 , f r , fnl , yl , aO , ratio , rb , tb ,btr , r2 , epp , f 1 , xl , zl) 
if (fl.gt.endl.and.fi. It.startl)  goto  25 
rmin=dsqr t ((a0-xl)*(a0-xl)+z l*z 1 ) 


Direction  of  normal  line  and  relative  velocity  at  pitch  point 


call  normal ( f 1 , fnl , xl , yl , zl , aO , ratio , rb , tb , btr , 
&  px,py,pz,rx,ry,rz) 

write(7,6)  fnl,xl,zl,px,py,pz 
write(8,6)  fnl,xl,zl,rx,ry,rz 
c 

nc=nc+l 
tp2(nc)=fn2 
fp2(nc)=fl 
pxj (nc)=xl 
pzj (nc)=zl 
xpn(nc)=px 
3q)n(nc)=py 
zpn(nc)=pz 
vx(nc)=rx 
vy(nc)=ry 
vz(nc)=rz 
c 
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25  do  30  i=l,ind 
xm(i)=0 
zm(i)=0 
rr(i)“0 
ff(i)-0 
30  continue 

C 

il  (iii 

c 

return 

end 


c 


c  Subroutine  WORM!  called  by  Subroutine  PITCH  and  ROOT 
c 


c  This  subroutine  is  used  to  calculate  the  coordinates  of  the  worm 
c  surface  in  Sj  specially  for  yj=0. 
c 

c  To  run  this  subroutine  the  input  parameters  are : 
c  aO  center  distance  of  the  first  enveloping  (mm) 

c  btr  apex  angle  of  the  generating  plane  (tad) 

c  ratio  speed  ratio  of  the  first  enveloping 

c  rb  radius  of  base  circle  (mm) 

c  tb  tan(bt) 

c  fnl  the  given  position  (tn2*ratio2)  (rad) 

c  ft  running  parameter  of  the  first  enveloping  (rad) 

c 
c 

c  The  output  of  the  subroutine  are: 

c  xm,zm  the  coordinates  of  the  worm  profile 

c  mw  the  index  of  the  worm  profile  point  on  edge  curve  (mw=l) 

c - - - - - - - - 

c 

subroutine  wormt (ft, fnl , y , aO , ratio , rb , tb , btr , xm , zm , mw) 
implicit  double  precision(a-h,o-z) 
c 

mw=0 

f2=ft/ratio 

fl=ft-fnl 

call  eqw(fl, f2 ,a0, ratio, rb,tb, 

&  p , q , r , s , pi , ql , rl , si , p2 , q2 , r2 , s2 ) 

c 

dw=p* (q2*rl - ql*r2 ) +q* (pl*r2 -p2*rl)+r* (p2*ql -pl*q2 ) 
if(dw.eq.O)  goto  10 

fz=- (p*(q2*sl-ql*s2)+q*(pl*s2-p2*sl)+s*(p2*ql-pl*q2) ) /dw 
dw=p*rl-pl*r 

xm=  (y*(ql*r-q*rl)+(sl*r-s*rl))/dw 
zm=  (y*(pl*q-p*ql)+(pl*s-p*sl))/dw 
if(zm.lt.fz)  goto  10 
goto  20 
c 

10  mw=l 
c 

20  return 
end 
c 

END  of  Subroutine  WORMT  ^ r - 

c 

c 
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c 

c 

c= 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


Subroutine  ROOT  called  by  Subroutine  PITCH 

This  subroutine  is  used  to  search  the  pitch  point  from  the  worm 
profile  in  Sj  specially  yj=0.  The  method  used  here  is  called 
golden  division. 


To  run  this  subroutine  the  input  parameters  are: 

aO  center  distance  of  the  first  enveloping  (mm) 

btr  apex  angle  of  the  generating  plane  (tad) 

r  pitch  radius  of  gear  (mm) 

ratio  speed  ratio  of  the  first  enveloping 
rb  radius  of  base  circle  (mm) 

tb  tan(bt) 

fnl  the  given  position  (tn2*ratio2)  (tad) 

fi,fe  interval  of  fl  for  root  searching  (tad) 

epp  convergence  criterion  for  root  searching  (tad) 


The  output  of  the  subroutine  are: 

xl,zi  the  coordinates  of  the  worm  pitch  point  in  Sj 
fl  the  value  of  running  parameter  at  pitch  point 


subroutine  root ( f i , f e , fnl , y , aO , ratio , rb , tb , btr , r , epp , f 1 , xl , zl ) 
implicit  double  precision(a-h,o-z) 


110 


120 

130 

150 

160 

170 


call  wormt(f i, fnl, y,a0, ratio, rb,tb, btr, xm,zm,mw) 

gi=ds  qrt((a0-xm)*( aO -xm) +zm*zm) - r 

call  wormt(fe,fnl,y,aO,ratio,rb,tb,btr,xm,zm,mw) 

ge=ds  qrt((a0-xm)*(a0-xm) +zm*zm) - r 

df=fe-fi 

if (df. It. epp)  goto  150 

ftl=fi+.382*df 

ft2=fi+.618*df 


call  wormt ( f tl , fnl , y , aO , ratio , rb , tb , btr , xm , zm , mw) 
gl=dsqrt( (a0-xm)*(a0-xm)+zm*zm) -r 

call  wormt(ft2, fnl, y,a0, ratio, rb,tb, btr, xm,zm,mw) 

g2=dsqrt(  (a0-xm)*(a0-xmHziii*zm)  -r 

if(gl*gi.gt.0)  goto  120 

fe=ftl 

ge=gl 

goto  110 

if(g2*gi.gt.0)  goto  130 

fi=ftl 

fe=ft2 


gi=gl 
ge=g2 
goto  110 
fi=ft2 


goto  110 
fl=(fi+fe)/2.0 

call  wormt ( f 1 , fnl , y , aO , ratio , rb , tb , btr , xl , zl , mw) 

goto  170 

mc=l 


return 

end 

—  END  of  Subroutine  ROOT 


127 


c  Subroutine  NORMAL  called  by  Subroutine  PITCH 

C 

C 


c  This  subroutine  is  used  to  calculate  the  propertities  of  the  worm 

t  liiitiiSi. 

4 

C 

c  To  run  this  subroutine  the  input  parameters  are : 
c  aO  center  distance  of  the  first  enveloping  (nun) 

c  btr  apex  angle  of  the  generating  plane  (rad) 

c  fl  the  value  of  running  parameter  at  pitch  point 

c  ratio  speed  ratio  of  the  first  enveloping 

c  rb  radius  of  base  circle  (mm) 

c  tb  tan(bt) 

c  xl,zl  the  coordinates  of  the  worm  pitch  point  in  Sj 

c 

c 

c  The  output  of  the  subroutine  are: 

c  wnx,wny,wnz  the  normal  direction  components  of 

c  pitch  point  in  Sj 

c  xw,yw,zw  the  relative  velocity  components  of 

c  pitch  point  in  Sj 

c . - . - - - - - 

c 

subroutine  normal ( f 1 , f nl , xl , yl , zl , aO , ratio , rb , tb , btr , 

&  wnx,wny,wnz,xw,yw,zw) 

implicit  double  precision(a-h, o-z) 
c 

f2=fl/ratio 
cf2=dcos(f2) 
sf2=dsin(f2) 
cf  l-=dcos  (  f  1  -  fnl ) 
s  f l“ds in ( f 1 - fnl ) 
c 

wnx=cfl*sf2+tb*sfl 

wny=tb*cfl-sfl*sf2 

wnz=-cf2 

rx=xl* ( 1 - aO/dsqrt (xl*xl+yl*yl) ) 
ry=y 1* ( 1 - aO/ds  qr t (xl*xl+y l*yl ) ) 
rz=zl 

xw=wn3r*rz  -wnz*ry 
yw=wnz*rx - wnx*rz 
zw=wnx*ry - wny*rx 
rlu=dsqrt  (xw^w+yw*yw) 
rld=dsqrt (xw*xw+yw*yw+zw*zw) 
rlead=dacos (rlu/rld) 
return 
end 
c 

- L— . end  of  Subroutine  ROOT 

c 

c 

0..---===::=— - - - ■!.: -  .... -  ■  - -  ■  - - 

c  Subroutine  LOAD  called  by  MAIN 
c 

c  This  subroutine  is  used  to  calculate  load  sharing  among  the  meshing 
c  teeth.  All  the  calculations  are  under  the  given  position  tn2 .  To 
c  drive  this  subroutine  the  subroutines  DIMEN,  CONTACT  and  PITCH  must 
c  be  called  first, 
c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 


c 

c 

c 


To 


run  this 
a2 
nc 

nzmax 

power 

r2 

ratio2 

rpm 

t2 

uf 

xl,zl 

pnx.pny, 

cont 


subroutine  the  input  parameters  are: 

center  distance  of  wormgear  (num) 

the  number  of  real  meshing  pairs 
number  of  meshing  teeth 

output  power  from  gear  shaft  (hw) 

pitch  radius  of  gear  (nm) 

speed  ratio  of  wormgear 

input  speed  from  worm  shaft  (tpni) 


the  position  of  each  meshing  pair 
friction  coefficient 
the  coordinates  of  pitch  points  in  Sj 
pnz  the  normal  direction  components  of 
pitch  point  in  Sj 

the  length  of  contact  lines  on  each  meshing  pair 


The  other  results  which  are  not  called  by  other  subroutines  are 
written  in  the  data  file  under  the  name  of  'load.dat' 


subroutine  load  ( a2,ratio2,r2, nzmax, power, rpm, uf,nc,t2, 
&  xl , zl , pnx , pny , pnz , cont , wtor , wtorf , f orn , f ornf ) 

implicit  double  precision(a-h,o-z) 

dimension  t2(nzmax) ,cont(nzmax) ,xl(nzmax) ,zl(nzmax) , 

&  pnx (nzmax) , pny (nzmax) , pnz (nzmax) , 

&  rei(lO) ,refi(10) ,ww(10) ,wwf(10) 

9  format(2x,6fl2.4) 


write (9,*)  ' 
write (9,*)  ' 
write (9,*)  ' 
&rque ' 

write(9,*)  ' 

e  t 

normal 

load  and  worm  torque 

distribution' 

f 

pair 

No.  1 

normal  load 

1  worm  1  to 

write (9,*)  ' 
&uf OO ' 
write (9,*)  ' 

& . ' 

pair 

No.  1 

uf=0  ufoO 

1  uf=0  1 

pi=4 . dO*datan(l . dO) 

grpm=rpm/ratio2 

torqeg=9 . 550*1 . 0e6*power/grpm 

calculation  of  the  sum  of  the  contact  lines 


dd=0 . 0 
ddf=0.0 
do  20  i=l,nc 

pn=ds qr t ( pnx ( i ) *pnx ( i ) +pny ( i ) *pny ( i ) +pnz ( i ) *pnz ( 1 ) ) 

xn=dabs (pnx ( i ) /pn) 

yn=dabs (pny ( i ) /pn) 

zn=dabs (pnz ( i ) /pn) 

rei(i)=zn*(a2-xl(i))-xn*zl(i) 

pv=dsqrt ( r2*r 2+ (xl ( i ) *ratio2 ) **2 ) 

ref i ( i) =uf *r2*r2/pv 

ww(i)=yn 

wwf ( i ) =yn+uf *xl ( i ) *ratio2/pv 
dd=dd+re i ( i ) *cont ( i ) 
ddf=ddf+(rei(i) -refi(i) )*cont(i) 

20  continue 


129 


c 

c  calculate  normal  load  and  worm  torque  distribution  (uf=0) 
c 

do  30  j=l,nc 

gfor=torqeg*cont(j )/dd 

wtor-gfor*ww(j)*xl(j) 

c 


c  with  friction 


c 

gforf=torqeg*cont(j ) /ddf 
wtorf=gforf*wwf (j )*xl(j ) 

write (9, 9)  t2(j )*ratio2 ,gfor,gforf ,wtor ,wtorf 
30  continue 
c 

return 

end 

c 

END  of  Subroutine  LDAD  -  ■  ■ 

c 
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Appendix  3  The  Principal  Relative  Curvatures 


Assume  that  a  surface  is  a  simple  surface.  At  a  point  Q  of  Z^  there 
exists  a  normal  n.  The  tangent  plane  of  the  surface  Z^  at  Q  is  perpendicular 
to  the  normal  n.  Any  plane  containing  the  normal  n  is  perpendicular  to  the 
tangent  plane.  This  plane  is  called  a  normal  plane  of  the  surface  Z^  at  point 
Q.  The  intersection  of  the  normal  plane  and  the  tangent  plane  is  a  straight 
line.  The  direction  of  the  straight  line  is  called  a  tangent  direction  of 
surface  Z^  at  point  Q.  Therefore,  at  a  point  on  Z^  there  is  only  one  normal 
direction  and  only  one  tangent  plane,  but  there  are  infinitely  many  tangent 
directions  and  infinitely  many  normal  planes. 


At  point  Q  the  intersection  of  surface  Z^  and  a  normal  plane  gives  a 
planar  curve,  correspondingly  the  intersection  of  the  normal  plane  and  the 
tangent  plane  gives  a  tangent  direction.  The  curvature  of  the  planar  curve  at 
point  Q  is  called  the  normal  curvature  of  surface  Zi  in  this  tangent 
direction.  Let  the  point  Q  be  given  by  the  position  vector  r(u,v),  where  u,  v 
are  two  independent  variables.  Then  the  normal  curvature  is: 

_  Ldu^  •*-2Mdudv+Ndv^  .  ,  v 

“  Edu^+2Fdudv+Gdv^  '■  '  ' 


where  E=ru^  ,  F=r„T^,  G=r^^,  L=r„„-n,  M=r„^*n,  N=r.(,^-n;  and  n  is  the  normal 
direction  given  by  r^xr^/l  r„xr^|  .  If  the  position  of  Q  is  given  by  r(x,y,z) 
where  z=z(x,y)  represents  the  equation  of  the  surface,  then  E,  F,  G,  L,  M,  N 
can  be  written  as : 


E  =  1  +  (-1^)^ 
dx 


F=  |£|£,  G=  1  +  (|^)^ 
ox  ay  ay 


Af  = 

DB  ‘ 


DB 


where 


(a. 2) 
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The  sign  of  the  normal  curvature  is  determined  by  the  normal  direction. 
If  the  normal  direction  points  to  the  center  of  curvature  center,  the  normal 

curvature  is  positive,  otherwise  it  is  negative. 


are  different.  The  principal  curvatures  are  the  extreme  values  of  the  normal 
curvatures.  The  corresponding  tangent  directions  are  called  the  principal 
directions  (Kreyszig,  1959). 


If  two  surfaces  and  22  meshing  at  a  contact  point  Q,  the 

normals  of  the  two  surfaces  at  the  contact  point  coincide.  Along  a  tangent 
direction  let  the  normal  curvatures  of  the  two  surfaces  be  and  then 

in  this  direction  the  relative  curvature  of  two  surfaces  is  defined  as; 

jr,  =  Jc»>-kr 

The  principal  relative  curvatures  are  the  extreme  values  of  the  relative 
curvatures ,  the  corresponding  directions  are  the  principal  relative 
directions.  If  the  two  surfaces  are  simple  surfaces,  there  are  two  principal 
relative  curvatures  and  two  principal  relative  directions. 


If  the  two  surfaces  contact  along  a  line,  the  direction  of  the  contact 
line  is  one  of  the  principal  relative  direction  and  the  corresponding 
principal  relative  curvature  is  zero.  The  other  principal  relative  direction 
is  perpendicular  to  the  direction  of  the  contact  line  and  the  principal 
relative  curvature  is  (Litvin,  1989;  Shen,  et  al.  1983): 

^  _  E2N2-2F2M2*G2L2  (a. 4) 

^  E^G^-Ff 


In  the  second  enveloping  process  the  gear  surface,  Eqs .  27,  is  generated 
by  the  worm  surface,  Eqs.  18.  In  this  process  the  two  surfaces  are  in  line 
contact.  Let  P=Ae^^^’F,  the  equations  of  the  gear  surface  can  be  written  as: 
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P  =  0 


D^P  =  0 
D.P  =  0 


(a. 5) 


By  solving  Eqs.  18,  the  relations  of  ZiCxj.yi.V'Cxi.yi) )  can  be  found.  A  point 
on  the  worm  surface  can  then  be  expressed  as  r^Cxi.yi.Zj^Cxj^.yi)) .  Further,  by 
solving  Eqs.  27,  the  relations  of  Z2(x2,y2),  V'(x2.y2).  ^(X2,y2)  can  be  found. 

A  point  on  the  gear  surface  can  be  written  as  r2(x2,y2,Z2(x2,y2) )  •  Then,  by 
using  Eq.  a. 4  the  principal  relative  curvature  is  obtained.  In  this  equation, 
one  should  make  the  following  substitutions: 

„  „.F^Fy^  (a. 6) 


(a. 7) 


Jp 

II 

+ 

•“l 

-  ^  ■  fi- 

*1 

•*1 

2  F  F 

■^*1  ^Zl 

F. 

*1 

_  ■F’x,*! 

.  W.=  - 

DA 

/  «i- 

DA 

/  ■**! 

DA 

p2  +  p2  P  P  P  ^  P  ^ 

F2 - - - ^2  p  2 


(a. 8) 


L2=  m  (  Pe^er  2l>c  P,^e^+ 

M2=  m  {  bd  Fiji -(be  +  cd  cr® 

N^=m  (dPPfiji-  2de  P^^^*  e^ 


(a. 9) 


where 


DA  =  4 


(a. 10) 
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dF 

p  .  S^F 

dg  ' 

^  dgdr 

(q. 

II 

yi.  Zi.  +1) 

D  “ 

dP 

■P«- 

ds  ' 

p  - 

ds  dt 

yi,  Zi< 

^2'  y2*  ^2' 

(a. 12) 


(a. 13) 
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